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FOREWORD 


This  report,  prepared  "by  staff  members  of  Arthur  D.  Little, 
Inc. , Acorn  Park , Cambridge,  Massachusetts,  is  the  final  techni- 
cal report  on  a design  study  to  determine  the  feasibility  of  a 
rotary-reciprocating  thermal  compressor  for  use  in  spaceborne, 
cryogenic  refrigeration  systems.  The  work  was  carried  out  under 
U.S.A.F.  Contract  F336l5-7^-C-3082  (Arthur  D.  Little,  Inc.,  Case 
No.  77056).  The  contract  was  in  support  of  Project  No.  6lU6, 
Task  6lU6  03  19,  and  was  made  possible  by  the  use  of  Laboratory 
Director's  Fronds.  The  work  was  administered  under  the  direction 
of  the  Air  Force  Flight  Dynamics  Laboratory,  Vehicle  Equipment 
Division,  with  Lt.  David  C.  Brubaker,  AFFDL/FEE,  as  Project 
Engineer. 


This  report  covers  work  from  April  to  August  1971*  and  was 
released  by  the  author  in  October  1971*  for  publication  as  a tech- 
nical report . 

This  technical  report  has  lean  reviewed  and  is  approved. 


ABSTRACT 


This  report  presents  the  results  of  a design  study  program  on  a 
rotary  reciprocating  compressor  which  uses  heat  as  the  source  of  driving 
energy.  The  objective  of  the  study  program  was  to  determine  the  feasi- 
bility of  using  a thermal  compressor  in  a spaceborne,  cryogenic  refrigera- 
tion system.  The  report  explains  the  principles  of  operation  of  a thermal 
compressor,  and  presents  the  equations  which  describe  its  behavior.  It 
also  describes  a computer  program  constructed  to  manipulate  these  equations, 
he  results  of  a number  of  runs  made  with  this  program  are  discussed,  and 
the  design  of  a thermal  compressor  to  meet  a specific  set  of  performance 
requirements  is  described.  Finally,  a comparison  is  made  between  two  types 
of  spaceborne,  closed-cycle,  cryogenic  refrigeration  systems— one  incorpo- 
rating a thermal  compressor  and  the  other  incorporating  an  electrically 
driven  compressor.  1 

It  is  concluded  that  a thermal  compressor  is  a viable  alternative  to  an 
electrically  driven  compressor  in  refrigeration  systems  where  heat  is  pre- 
fer5!dj°VeL  electricity  as  the  energy  source.  For  the  specific  case 
studied,  the  refrigeration  system  incorporating  a thermal  compressor  required 
92  watts  of  electrical  power  and  1,064  watts  of  heat  at  1,250°F;  while  the 
system  incorporating  an  electric  compressor  required  587  watts  of  electri- 
cal  power.  The  weights  and  sizes  of  the  two  refrigeration  systems  were 

comparable.  No  serious  technical  problems  are  foreseen  in  building  a 
thermal  compressor. 
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I.  INTRODUCTION  AND  SUMMARY 


I For  a number  of  years,  the  Air  Force  has  supported  the  development 

of  rotary-reciprocating  refrigeration  systems.  These  systems  have  been 
designed  to  produce  refrigeration  at  cryogenic  temperatures  in  spaceborne 
applications  where  long,  maintenance— free  life  and  low  input  power  are 
primary  system  requirements.  A complete  77  K refrigeration  system  has  been 
developed  (see  Reference  1)  and  components  of  a 3.6  K refrigerator  have 
been  built  and  tested  (see  Reference  2).  A system  which  produces  refrigera- 
tion simultaneously  at  12  K and  60  K is  presently  under  development  (see 
Reference  3) . All  of  these  systems  utilize  the  reversed  Brayton  thermo- 
dynamic cycle  and  employ  reciprocating  compressors  and  expansion  engines. 

In  all  these  systems,  all  of  the  input  power  to  the  refrigerator  has  been 
in  the  form  of  electrical  energy — with  approximately  80Z  of  this  energy 
being  used  to  compress  the  gas  in  the  compressor  working  volumes,  and  the 
balance  being  used  to  rotate  pistons,  to  overcome  losses,  and  to  actuate 
the  expander  pistons. 

In  recent  years,  the  Air  Force  has  become  interested  in  reducing  the 
electrical  power  requirements  of  spaceborne,  closed-cycle  refrigeration 
systems  by  using  heat  rather  than  electrical  energy  as  the  major  energy 
source  for  the  refrigerator.  In  a refrigerator  which  uses  the  reversed 
Brayton  cycle,  a major  reduction  in  the  electrical  power  requirement  can 
be  accomplished  by  employing  a compressor  which  uses  heat  rather  than 
electricity  to  supply  the  work  of  compression.  A heat-driven  compressor 
of  this  sort,  generally  referred  to  as  a thermal  compressor,  may  be  used 
in  place  of  the  electrically  driven  compressor  with  relatively  minor 
changes  to  other  portions  of  the  refrigeration  system. 

A thermal  compressor  is  a reciprocating  machine  which  consists  of  two 
basic  elements:  a cylinder  assembly  and  a drive  mechanism.  The  cylinder 

assembly  consists  of  a cylinder  and  a displacer.  The  displacer  separates 
the  cylinder  into  two  working  volumes.  These  two  working  volumes  are 
connected  by  flow  passages  which  contain  heat  exchangers  and  a thermal 
regenerator.  One  end  of  the  cylinder  is  maintained  at  an  elevated  tempera- 
ture by  the  addition  of  heat,  and  the  other  end  is  maintained  at  ambient 
temperature  by  rejecting  heat  to  a heat  sink.  As  the  displacer  is  moved 
back  and  forth  in  the  cylinder  by  the  drive  mechanism,  the  pressure  in  the 
cylinder  fluctuates.  Valving  in  the  ambient  temperature  working  volume 
communicates  this  working  volume  with  inlet  and  discharge  manifolds  when 
the  pressure  in  the  cylinder  matches  the  pressure  in  the  manifolds.  All 
of  these  actions  result  in  gas  being  compressed  from  inlet  pressure  to 
discharge  pressure  by  the  heat  that  is  added  to  the  hot  end  of  the  cylinder. 
The  principles  of  operation  of  a thermal  compressor  are  described  in  more 
detail  in  Section  II  of  this  report. 


In  a rotary-reciprocating  thermal  compressor,  the  displacer  is  actuated 
by  an  electrically-powered  drive  assembly  which  is  similar  to  that  used  in 
an  electric  rotary-reciprocating  compressor.  The  drive  assembly  comprises 
a rotary  motor  to  rotate  the  displacer,  a linear  actuator  to  reciprocate  the 
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displacer,  and  a gas  spring  to  provide  a restoring  force  to  the  displacer. 

The  rotation  activates  self-acting  gas  bearings  which  completely  support 
the  displacer.  Gas  flows  In  and  out  of  the  compressor  are  regulated  by 
ports  In  the  cylinder  assembly  which  are  opened  and  closed  by  the  rotating, 
reciprocating  action  of  the  displacer.  The  electrical  power  required  by  the 
drive  assembly  Is  determined  prlnlcpally  by  the  power  required  to  move  the 
displacer.  Since  the  pressures  at  both  ends  of  the  displacer  are  approximately 
the  same,  relatively  little  electrical  energy  Is  required  to  move  the  displacer. 
Thus,  the  electric  Input  to  the  compressor  Is  minimal,  the  major  source  of 
Input  energy  being  the  heat  supplied  to  the  hot  end  of  the  cylinder. 

In  recognition  of  the  potential  advantages  of  a thermal  compressor  In 
reducing  electrical  Input  power  to  spacebome  refrigeration  systems,  the 
Air  Force  has  supported  the  design  study  program  which  is  summarised  In 
this  report.  The  objective  of  the  study  program  has  been  to  determine  the 
feasibility  of  a thermal  compressor  for  a rotary-reciprocating  cryogenic 
refrigeration  system.  The  study  has  consisted  of  analyzing  and  designing 
a thermal  compressor  to  the  point  where  technical  feasibility  and  characteris- 
tics could  be  determined,  and  then  comparing  a rotary-reciprocating  refri- 
gerator Incorporating  a thermal  compressor  with  a rotary-reciprocating 
refrigerator  Incorporating  an  electrical  compressor — with  both  systems  being 
designed  to  meet  a specified  set  of  refrigeration  requirements.  The  refri- 
geration requirements,  which  were  specified  by  the  Air  Force  at  the  outset 
of  the  study,  are  listed  In  Table  1. 


TABLE  I 


REFRIGERATION  REQUIREMENTS 


Refrigeration  Capacity 

0.3  Watts  at  12  K 

3.0  Watts  at  28  K 

6.0  Watts  at  75  K 
Heat  Rejection  Temperature 


300  K 


The  first  step  In  the  study  program  was  to  determine  an  appropriate 
thermodynamic  cycle  for  a refrigeration  system  employing  a thermal  compressor 
This  was  done  to  determine  the  required  performance  specifications  for  the 
thermal  compressor  and  to  Identify  ranges  of  cycle  parameters  that  would 
lead  to  practical  designs  of  the  compressor.  This  task  was  Included  in 
recognition  of  the  fact  that  the  most  suitable  thermodynamic  cycle  utilizing 
a thermal  compressor  might  have  a different  pressure  ratio  and  pressure 
level  than  If  the  cycle  Incorporated  an  electrically  driven  compressor. 

The  cycle  selection  process  is  discussed  In  Section  III  of  this  report,  and 
the  compressor  specifications  resulting  from  the  cycle  selected  are  listed 
In  Table  II. 


TABLE  II 


SPECIFICATIONS  FOR  THE  THERMAL  COMPRESSOR 


Process  Fluid 
Mass  Flow  Rate 
Inlet  Pressure 
Discharge  Pressure 
Inlet  Temperature 


Helium  Gas 
3.66  LB/HR 
29.4  LB /IN 2 
58.8  LB/IN2 
300  K 


A parametric  analysis  was  conducted  in  order  to  design  the  working 
volumes  of  the  thermal  compressor  and  to  determine,  in  a general  way, 
the  characteristics  of  a thermal  compressor.  The  analysis  consisted  of 
deriving  the  equations  which  describe  the  various  elements  of  the  com- 
pressor (working  volumes,  heat  exchangers,  regenerators,  etc.).  These 
equations  were  then  programmed  for  a digital  computer  and  a number  of 
cases  were  run  with  the  program.  The  basis  for  the  parametric  analysis 
and  the  results  of  the  analysis  are  presented  in  Section  IV  of  this 
report.  A listing  of  the  equations  which  describe  a thermal  compressor 
is  contained  in  Appendix  I and  the  computer  program  based  on  these 
equations  is  described  in  Appendix  II. 


I Based  on  the  results  of  the  parametric  a.ialysis,  we  developed  the 

preliminary  design  of  a thermal  compressor  to  meet  the  performance  speci- 
fications which  were  listed  in  Table  II.  This  design  is  summarized  in 
Table  III,  and  an  outline  drawing  of  the  machine  is  shown  in  Figure  1. 

The  design  of  the  thermal  compressor  is  discussed  in  more  detail  in  Section  V. 


■ 


An  overall  comparison  between  a rotary-reciprocating  refrigerator 
employing  a thermal  compressor  and  one  employing  an  electric  compressor  is 
shown  in  Table  IV.  From  this  table  it  can  be  seen  that  a refrigerator  uti- 
lizing a thermal  compressor  is  somewhat  heavier  than  a refrigerator  incor- 
porating an  electric  compressor.  The  refrigerator  incorporating  a thermal 
compressor  requires  495  watts  less  electrical  power  than  an  electric-driven 
refrigerator,  but  will  require  1,064  watts  of  thermal  power  at  1,250°F. 

The  two  refrigerators  are  virtually  the  same  size.  A more  detailed  compari- 
son between  the  two  refrigeration  systems  is  contained  in  Section  VI. 

We  conclude  that  the  thermal  compressor  is  a viable  alternative  to  an 
electrically  driven  compressor  in  spaceborne,  rotary  reciprocating,  cryogenic 
systems,  and  recommend  that  thermal  compressor  refrigerators  be  considered 
along  with  electric  compressor  refrigerators  in  spacecraft  system  studies. 
Further  conclusions  drawn  from  the  study  are  contained  in  Section  VII. 
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TABLE  III 


SUMMARY  OF  THERMAL  COMPRESSOR  DESIGN 


3.66  LB/HR 

29.4  LB/IN^  (2  Atmospheres) 
58.8  LB/IN^  (4  Atmospheres) 
300  K 


Inlet  Pressure 


Discharge  Pressure 


Inlet  Temperature 

Number  of  Stages 

Number  of  Cylinders  Per  Stage 

Compressor  Size 


See  Figure  1 


Compressor  Weight 


Input  Power 


Thermal  (At  1,250°F) 
Electrical  (100  V DC) 


1(064  Watts 


The  system  may  also  be  designed  for  28  V DC 


TABLE  IV 


COMPARISON  OF  REFRIGERATOR  TYPES 


Refriserator 

Refriaerator 

with  Thermal 

with  Electric 

Compressor 

Compressor 

WEIGHT 

Compressor  Assembly 

155  LB 

70  LB 

Expander  Assembly 

70 

70 

Power  Conditioning  Equipment 

30 

35 

TOTAL 

255  LB 

175  LB 

INPUT  POWER 

Thermal  (At  1,250°F) 

1,064  Watts 

— 

Electrical  (100  V DC) 

92  Watts 

587  Watts 

II.  PRINCIPLES  OF  OPERATION  OF  A THERMAL  COMPRESSOR 


A thermal  compressor  comprises  two  basic  subelements:  a cylinder 

which  contains  a reciprocating  displacer,  and  a drive  assembly  which  Im- 
parts motion  to  the  displacer.  The  cylinder  consists  of  a single,  fixed 
volume  which  Is  divided  Into  two  working  volumes  by  the  displacer  (see 
Figure  2) . One  working  volume  Is  maintained  at  a high  temperature  (typi- 
cally approximately  1,000°F)  by  the  addition  of  heat,  and  the  other  working 
volume  Is  maintained  at  ambient  temperature  by  the  removal  of  heat.  The 
two  working  volumes  communicate  with  each  other  through  a thermal  regen- 
erator. The  total  internal  volume  of  the  device  Includes  void  volumes 
in  addition  to  the  two  working  volumes.  The  regenerator  may  be  incorpo- 
rated either  in  the  displacer  (as  in  Figure  2)  or  it  may  be  mounted  in  the 
cylinder.  Insulation  surrounds  the  hot  end  of  the  device,  and  the  design 
is  such  that  heat  losses  from  the  hot  end  are  minimized.  The  ambient 
working  volume,  is  connected  to  two  manifolds:  a sucMon  (inlet)  manifold 

and  a discharge  manifold.  The  connection  is  through  valving  or  through 
porting,  and  is  indicated  schematically  in  Figure  2 by  check  valves. 

As  the  displacer  moves  back  and  forth,  the  working  volumes  undergo 
pressure-volume  events  characteristic  of  reciprocating  compressors  and 
expansion  engines.  The  resulting  idealized  pressure-volume  (P-V)  diagrams 
are  shown  in  Figure  2.  The  P-V  diagram  of  the  hot  working  volume  is  char- 
acteristic of  an  isothermal  expansion  engine,  with  the  work  output  per 
cycle  being  equal  to  the  heat  added  per  cycle . This  work  is  delivered  to 
the  displacer,  and  thence  to  the  gas  in  the  ambient  working  volume.  The 
ambient  working  volume  has  a P-V  diagram  characteristic  of  an  isothermal 
compressor.  In  this  P-V  diagram,  the  work  input  to  the  gas  per  cycle  is 
equal  to  the  heat  rejected  from  the  cylinder  each  cycle.  Since  the 
expander  work  and  the  compressor  work  are  equal  to  each  other,  the  heat 
rejection  at  the  ambient  end  is  equal  to  the  heat  added  at  the  hot  end. 

Gas  is  taken  from  the  inlet  manifold  into  the  compressor  during  the 
process,  b-c,  and  this  gas  is  delivered  to  the  higher  pressure  discharge 
manifold  during  the  process,  d-a.  In  the  idealized  case  (no  pressure  drop 
and  no  friction)  no  mechanical  work  is  required  to  move  the  displacer  since 
the  pressure  is  the  same  at  each  end  of  the  displacer.  Therefore,  the  only 
form  of  input  energy  is  the  heat  added  to  the  hot  end  of  the  device.  Thus, 
the  compressor  utilizes  heat  to  compress  gas  from  suction  pressure  to  dis- 
charge pressure.  It  is,  in  effect,  a heat  engine  (operating  between  the 
hot  temperature  and  the  ambient  temperature)  driving  an  isothermal  compressor 
operating  at  the  ambient  temperature. 

In  a practical  device,  there  are  pressure  drop  losses.  The  drive 
assembly  must  supply  work  to  the  displacer  to  overcome  these  losses.  As 
will  be  seen  in  a following  section,  the  presence  of  the  drive  shaft  which 
connects  the  ambient  end  of  the  displacer  to  the  drive  assembly  slightly  re- 
duces the  amount  of  work  delivered  to  the  compressor  working  volume  from 
the  expander  working  volume,  since  the  areas  of  the  two  working  volumes 
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of  Discharge 


are  different.  This  results  in  the  expander  working  volume  delivering 
some  work  to  the  displacer  assembly,  and  thus  reduces  the  amount  of  work 
which  must  be  provided  by  the  drive  assembly. 

In  a rotary-reciprocating  thermal  compressor,  the  drive  assembly 
rotates  the  displacer  as  well  as  reciprocating  it.  The  rotation  activates 
self-acting  gas  bearings  which  support  the  displacer.  Further,  ports  in 
the  drive  shaft  are  opened  and  closed  by  the  displacer  motion,  to  control 
gas  flows  into  and  out  of  the  ambient  working  volume. 

The  thermal  compressor  is  not  a new  idea,  nor  does  it  require  develop- 
ment of  new  technology.  A thermal  compressor  was  patented  by  Bush  in 
1939  (Reference  4) . In  recent  years , the  National  Institutes  of  Health 
have  sponsored  several  developments  of  miniature  thermal  compressors  for 
use  as  a component  in  a completely  implantable  artificial  heart  system 
(Reference  5,  6,  and  7).  Technology  directly  applicable  to  thermal  compres- 
sors have  been  sponsored  by  the  U.S.  Air  Force  and  NASA  in  the  form  of 
development  of  Vuilleumier  (VM)  cycle  refrigerators  (the  hot  end  of  the 

refrigerator  Is,  in  effect,  a thermal  compressor  which  supplies  compressed 
gas  to  the  cold  and  of  the  refrigerator).  Additionally,  the  technology 
associated  with  Stirling  cycle  engines  is  also  applicable  to  the  thermal 
compressor.  Thus,  there  is  a technology  base  on  which  to  draw  for  the 
development  of  a rotary  reciprocating  thermal  compressor. 
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III.  THERMODYNAMIC  CYCLE  ANALYSIS 


The  thermodynamic  cycle  analysis  has  been  performed  by  using  an 
existing  rotary  reciprocating  computer  program  to  analyze  a number  of 
reversed  Bray ton  cycles  (see  Reference  8) . The  object  of  the  analysis 
has  been  to  select  a thermodynamic  cycle  for  a refrigerator  incorporating 
a thermal  compressor,  and  to  set  the  performance  of  specifications  for 
the  compressor. 

Our  approach  to  the  cycle  selection  has  been  to  use  the  system 
analysis  computer  program  to  determine  the  characteristics  of  an  electri- 
cally driven  refrigerator  as  a function  of  various  cycle  parameters.  For 
purposes  of  sizing  a thermal  compressor,  we  have  selected  the  cycle  which 
has  the  "best"  combination  of: 

1.  Low  input  power  to  the  electrically  driven  compressor 

2.  Low  system  weight  for  an  electrically  driven  compressor 

3.  Reasonable  pressure  ratio  for  the  thermal  compressor 

4.  Reasonable  expander  assembly  design 

The  basis  for  using  the  results  of  cycle  analysis  on  the  electrically 
driven  compressor  for  selecting  a cycle  for  the  thermal  compressor  is  that 
the  trend  of  input  power  to  the  thermal  compressor  will  follow  the  trend 
of  input  power  to  the  electrically  driven  compressor.  Thus,  the  cycle 
which  results  in  minimum  power  for  the  electrically  driven  compressor  will 
also  result  in  minimum  input  power  to  the  thermal  compressor.  Since  re- 
frigerator weight  is  a direct  function  of  input  power,  this  selection  process 
will  also  result  in  a refrigerator  of  minimum  weight. 

In  the  cycle  analysis,  most  cycle  parameters  can  be  considered  to  be 
constant,  with  only  a few  carried  as  variables.  The  fact  that  there  are 
three  heat  loads  dictates  a three  expander  cycle  incorporating  five  process 
heat  exchangers.  The  following  values  have  been  used  to  describe  the  ex- 
panders and  heat  exchangers: 

Efficiency  of  Expander  No.  1 80Z 

Efficiency  of  Expander  No.  2 75Z 

Efficiency  of  Expander  No.  3 70Z 

Heat  exchanger  effectiveness  98Z  (all  exchangers) 


Heat  exchanger  pressure  loss 
coefficient 


.01  (all  exchangers) 


These  values  are  representative  of  values  obtainable  In  practice. 
There  are  only  two  variables  In  addition  to  those  mentioned  above 
required  In  order  to  completely  specify  the  thermodynamic  cycle.  They 
are  the  two  pressure  levels  In  the  cycle.  A number  of  cases  have  he,en 
run,  considering  the  two  presuure  levels  as  variables.  The  cases  which 
have  been  run  are  within  the  following  range: 


Low  side  pressure 
Pressure  ratio ^ 


.5  to  4 atmospheres 


1.5  to  4 


A refrigerator  design  could  not  be  obtained  for  pressure  ratios  less  than 
2.0  because  the  work  extraction  capability  of  the  expanders  was  not  suf- 
ficient to  offset  the  parasitic  heat  loads.  As  a result,  the  cycle  selec- 
tion Is  based  on  pressure  ratios  between  2 and  4. 

The  refrigerator  characteristics  of  primary  Interest  to  the  cycle 
selection  process  are  the  refrigerator  Input  power  and  the  refrigerator 
weight.  These  two  characteristics  are  plotted  In  Figures  3 and  4 as  a 
function  of  cycle  pressure  ratio  and  cvcle  low  side  pressure.  It  should 
be  noted  that  the  weights  shown  In  Figure  4 are  lower  than  may  actually 
be  achievable  In  practice,  but  the  trends  are  correct. 

Based  on  these  curves,  we  draw  the  following  conclusions: 

1.  The  minimum  Input  power  Is  obtained  at  a cycle  pressure 
ratio  of  2.5. 

2.  The  pressure  ratio  at  which  the  minimum  Input  power  is 
obtained  is  relatively  Insensitive  to  the  low  side  pressure 
(i.e.,  the  cycle  pressure  level). 

3.  The  refrigerator  weight  decreases  with  Increasing  cycle 
pressure  level  and  with  Increasing  cycle  pressure  ratio. 

Note  on  Figure  3 that  the  vertical  axis  is  compressor  shaft  power. 

This  Is  the  shaft  power  that  would  have  to  be  supplied  to  an  elec- 
trically driven  compressor  which  has  an  lsentroplc  efficiency  of  84%. 

This  power  figure  takes  into  account  Inefficiencies  In  the  compression 
processes,  but  not  losses  (e.g.,  in  the  gas  springs,  linear  actuators, 
power  supplies,  etc.),  or  power  to  drive  the  rotary  motors. 

Based  solely  on  these  conclusions,  we  would  select  a thermodynamic 
cycle  with  a pressure  ratio  of  2.5  (to  minimize  Input  power)  and  a low 


Pressure  Ratio  - 


high  side  (discharge)  pressure 
low  side  (suction)  pressure 


320 


84%  (Uentropic) 
Compressor  Efficiency 


side  pressure  level  of  3-4  atmospheres  (to  minimise  the  weight  of  the 
refrigerator).  However,  the  practical  aspects  of  refrigerator  design  must 
also  be  considered  in  selecting  a thermodynamic  cycle.  The  low  temperature 
(12  K)  expansion  engine  has  a very  low  capacity,  and  will  have  a small  diameter 
piston  in  any  event.  This  engine  can  be  made  larger  if  we  reduce  the  pres- 
sure ratio  and  pressure  level  of  the  cycle.  Since  a pressure  ratio  of  2, 
and  a low  side  pressure  of  2 atmospheres  results  in  a reasonably  sized 
(i.e.,  .5  inch  diameter  piston)  low  temperature  expansion  engine,  these 
parameters  were  selected  for  the  thermodynamic  cycle.  The  penalties  of 
selecting  this  cycle  are  not  severe.  In  reducing  the  pressure  ratio  to 
2 from  the  optimum  pressure  ratio  of  2.5,  we  increase  the  compressor 
shaft  power  only  slightly  (see  Figure  3).  Additionally,  we  do  not  change 
the  size  of  the  displacer  in  the  thermal  compressor,  since  the  increase 
in  gas  flow  rate  required  by  the  lower  pressure  ratio  cycle  is  offset  by 
the  increase  in  volumetric  efficiency  of  the  compressor  at  the  lowev 
pressure  ratio  (see  Section  IV,  B) . 

Selection  of  a three  expander  cycle  is  an  obvious  choice  from  a 
thermodynamic  standpoint  because  work  is  extracted  from  the  gas  at  tempera- 
tures very  near  those  at  which  the  heat  loads  are  absorbed.  An  alternative 
to  the  three  expander  cycle  is  a two  expander  cycle  in  which  the  second 
stage  engine  absorbs  both  the  lower  temperature  heat  loads  at  12  K.  The 
rationale  for  considering  a cycle  of  lower  thermal  efficiency  is,  first  of 
all,  that  two  expanders  rather  than  three  would  be  required,  and  secondly 
that  the  second  stage  engine  would  be  larger  (the  bore  would  be  about  0.75 
inch)  and  therefore  more  easily  designed. 

We  analyzed  a two  expander  cycle  with  pressure  ratios  between  1.5  and 
4.0  and  with  low  side  pressures  between  0.5  atmospheres  and  4.0  atmospheres. 
Qualitatively,  the  results  were  the  same  as  in  the  three  expander  cycle, 
viz.,  the  minimum  input  power  is  obtained  at  a pressure  ratio  of  2.5  and 
is  relatively  insensitive  to  low  side  pressure.  However  the  compressor 
shaft  power  for  the  two  cycle  types  were  vastly  different — the  input  power 
for  a refrigerator  utilising  the  two  expander  cycle  being  over  three  times 
that  for  a refrigerator  utilizing  the  three  expander  cycle.  With  a com- 
pressor efficiency  of  84%,  the  two  expander  cycle  compressor  would  require 
929  watts  of  input  power  compared  to  307  watts  for  a three  expander  cycle 
compressor.  This  increase  in  input  power  far  outweighs  any  advantages 
which  might  be  gained  in  the  mechanical  design  of  the  two  expander  refri- 
gerator. Consequently,  the  two  expander  cycle  was  not  considered  further. 


The  thermodynamic  cycle  selected  for  the  thermal  compressor  refrigerator 
is  shown  in  Figure  5.  The  thermal  compressor  discussed  in  following  sections 
has  been  designed  to  meet  the  requirements  of  this  cycle. 


IV.  PARAMETRIC  ANALYSIS  OF  THERMODYNAMIC  PROCESSES 


A.  GENERAL 

This  section  discusses  the  parametric  analysis  of  the  thermodynamic 
processes  which  occur  in  a thermal  compressor,  and  also  presents  the 
results  of  the  analysis  and  the  conclusions  drawn  from  it.  We  present 
a physical  description  of  the  analytical  model  used  to  define  the  ther- 
mal  compressor,  and  the  assumptions  upon  which  the  model  is  based.  In 
the  course  of  the  analysis,  we  derived  a number  of  equations  which  des- 
cribe the  thermodynamic  processes  occurring  in  the  compressor  and  which 
are  used  in  arriving  at  a design  and  determining  its  characteristics. 
These  equations  are  listed  and  discussed  in  Appendix  I.  Since  a large 
number  of  equations  is  required  to  define  a thermal  compressor,  and 
since  there  are  a number  of  independent  variables  to  consider  in  the 
design  optimization  process,  we  constructed  a digital  computer  program 
to  manipulate  the  equations.  This  program  is  covered  briefly  in  this 
section,  and  is  described  in  detail  in  Appendix  II.  All  terminology 
used  in  this  section,  as  well  as  in  Appendices  I and  II  is  defined  in 
the  list  of  abbreviations  and  symbols.  The  computer  program  was  used 
extensively  in  the  course  of  designing  a thermal  compressor  to  meet  the 
performance  specifications  which  resulted  from  the  thermodynamic  cycle 
analysis.  The  final  part  of  this  section  presents  thn  results  of  these 
computer  runs,  and  the  conclusions  drawn  from  the  results.  The  section 
concludes  with  a description  of  a baseline  design  point  to  which  the 
thermal  compressor  was  designed. 


ANALYTICAL  MODEL 
1.  Operating  Characteristics 


This  analysis  is  concerned  with  the  processes  occurring  in  the 
hot  and  ambient  ends  of  the  cylinder,  in  the  regenerator,  in  the  hot 
and  ambient  heat  exchangers  and  in  void  volumes.  It  is  also  concerned 
with  how  these  processes  interact  to  determine  the  pressure  ratio/ flow 
performance  characteristics  and  the  power  input  requirements  to  the 
cylinder  and  the  displacer. 

It  is  convenient  to  subdivide  the  cyclic  operation  of  the  compres- 
sor into  four  regimes:  re-expansion,  suction,  compression,  and  dis- 

charge. The  positions  of  the  displacer  at  the  four  points  corresponding 
to  the  transition  from  one  regime  to  the  next  are  shown  in  Figure  6.  A 
diagram  of  the  pressures  in  the  hot  and  ambient  ends  versus  stroke,  on 
which  the  four  -positions  are  indicated  as  _a,  b^  c and  <1,  is  shown  in 
Figure  7.  At  Position  <*,  most  of  the  gas  in  the  cylinder  is  at  the  hot 
end.  During  the  re-expansion  process  from  £ to  l>,  both  the  inlet  and 
discharge  ports  in  the  ambient  end  are  closed,  isolating  the  cylinder 
from  the  suction  and  discharge  reservoirs.  Tke  displacer  moves  toward 
the  hot  end,  displacing  gas  from  there  to  the  ambient  end  and  causing 
the  pressure  in  the  cylinder  to  decrease  from  the  discharge  pressure 
level  to  the  suction  pressure  level.  At  Point  b,  when  the  suction 
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Hot  Haat  Exchanger 


Ambiant  Haat  Exchanger 


a.  End  of  Discharge 


b.  Start  of  Suction 


c.  End  of  Suction 


d.  Start  of  Discharga 


FIGURE  6 DISPLACER  POSITIONS 


pressure  Is  reached,  the  suction  port  opens,  admitting  gas  to  the 
cylinder  from  the  suction  line.  Movement  of  the  displacer  continues 
from  b to  c until  the  displacer  reaches  Its  extreme  position  near  the 
hot  end,  at  which  time  the  suction  stroke  Is  complete.  In  the  compres- 
sion process,  from  c_  to  <1,  again  both  the  suction  and  discharge  ports 
are  closed,  and  gas  is  displaced  from  the  ambient  end  to  the  hot  end, 
causing  the  pressure  in  the  cylinder  to  rise.  At  Point  <1  the  discharge 
pressure  level  is  reached,  the  discharge  port  opens  and  flow  from  the 
ambient  end  of  the  cylinder  through  the  discharge  line  begins.  During 
the  discharge  process,  <1  to  _a,  the  pressure  level  in  the  cylinder  is 
maintained  constant  by  the  combination  of  mass  discharge  from  the 
ambient  end  and  gas  displacement  from  the  ambient  end  to  the  hot  end. 

The  pressures  in  both  the  hot  and  ambient  ends  are  shown  in  Fig- 
ure 7.  The  pressure  difference  between  them  is  caused  by  pressure  drop 
in  the  regenerator  and  in  the  ambient  and  hot  heat  exchangers.  Its 
direction  is  determined  by  the  direction  of  flow  through  the  regenera- 
tor. From  Figure  7,  one  can  see  that  the  pressure  versus  stroke  diagram 
of  the  ambient  end  is  similar  to  that  of  a compressor,  while  the  diagram 
of  the  hot  end  resembles  that  of  an  expander. 

Volume  definitions  and  sign  conventions  used  in  the  analysis  are 
shown  in  Figure  8.  The  displacement  volumes  at  the  hot  and  ambient  ends 
are  somewhat  different  because  of  the  shaft  connected  to  the  ambient  end 
of  the  displacer.  The  void  volumes  indicated  in  Figure  8 (l.e.,  V2V, 
V3V,  and  V4V)  are  comprised  of  the  various  elements  summarized  in 
Table  V. 


2.  Cylinder  Processes 

In  analyzing  the  processes  occurring  in  the  cylinder,  we  have  made 
the  following  assumptions: 

a.  Isothermal  conditions  exist  in  the  gas  in  the  hot  and 
ambient  ends. 

b.  Inlet  and  discharge  porting  pressure  drops  are  negligible 
near  the  ends  of  the  displacer  stroke. 

c.  Temperature  of  gas  in  void  volumes  will  be  either  TH,  TC 
or  an  intermediate  temperature  equal  to  the  log  mean 
temperature. 

d.  The  displacer  motion  is  sinusoidal. 

e.  The  working  fluid  is  a perfect  gas. 

Assumption  1 is  usually  made  in  analyses  of  this  class  of  device, 
i.e.,  in  Vuilleumier  cycle  and  Stirling  cycle  analyses.  During  the  con- 
stant pressure  suction  and  discharge  processes,  the  temperature  of  the 
gas  in  the  hot  and  ambient  ends  should,  in  fact,  be  isothermal.  During 
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FIGURE  8 VOLUME  DEFINITIONS 


TABLE  V 


VOID  VOLUME  SUMMARY 


HOT  END  (V3V) 

* Fluid  Passage,  Cylinder  to  Hot  Heat  Exchanger 

* Hot  Heat  Exchanger 

* Cylinder  End  Clearance 


AMBIENT  END  (V2V) 

• Fluid  Passage t Cylinder  to  Ambient  Heat  Exchanger 

• Ambient  Heat  Exchanger 

• Cylinder  End  Clearance 

• Inlet /Discharge  Porting  Passages 


INTERMEDIATE  (V4V) 

* Regenerator 

* Displacer/Cylinder  Radial  Clearance 


the  re-expanslon  process,  there  Is  a tendency  for  the  gas  In  the 
cylinder  to  cool;  and  during  the  compression  process,  there  Is  a ten- 
dency for  the  gas  to  warm  up.  Maintaining  Isothermal  conditions  during 
the  latter  two  processes  will  rely  on  adequate  heat  transfer  to  and 
from  the  cylinder  walls.  Such  heat  transfe'  will  be  augmented  by  cir- 
culations In  the  gas  Induced  by  rotation  of  the  piston.  For  devices  In 
the  size  range  of  our  Interest  and  at  the  low  operating  speeds  that  are 
appropriate  for  a thermal  compressor,  we  believe  that  Assumption  1 Is  a 
reasonable  approximation. 

Assumption  2 Is  based  on  the  fact  that,  with  sinusoidal  displacer 
motion,  the  flow  rate  Into  the  cylinder  as  Point  c_  Is  approached,  and 
out  of  the  cylinder  as  Point  a_  Is  approached,  goes  to  zero.  Hence, 
porting  pressure  drop  should  be  negligible  at  the  ends  of  the  stroke. 
Experience  with  the  electric-driven  rotary-reciprocating  compressor, 
Including  measurements  of  pressure-volume  diagrams,  has  Indicated  that 
this  assumption  Is  valid. 

Assumptions  3 and  5 are  good  approximations  and  are  self-explanatory. 
Assumption  4 Is  also  an  excellent  approximation  since  the  oscillating 
displacer  mass  experiences  only  a small  damping  force  due  to  the  pres- 
sure difference  between  the  hot  and  ambient  ends. 

There  Is  a pressure  difference  between  hot  and  ambient  ends  which 
Is  caused  by  pressure  drop  through  the  regenerator  and  the  heat  ex- 
changers. It  determines  the  necessary  electric  power  Input  to  maintain 
motion  of  the  displacer.  Our  approach  has  been  to  calculate  the  pres- 
sure difference  between  hot  and  ambient  ends  as  a function  of  the  dis- 
placer position  and  to  Integrate  the  resulting  force  times  differential 
displacement  throughout  a cycle  to  obtain  the  displacer  power  Input. 

Using  the  assumptions  above,  It  Is  possible  to  determine  the 
relationships  between  cylinder  pressure  and  the  mass  of  gas  In  each  of 
the  volumes  as  a function  of  displacer  position.  These  relationships 
lead  to  definition  of  cylinder  size,  power  input  and  other  Important 
parameters  required  to  achieve  given  flow  rate  and  suction  and  discharge 
pressures.  The  equations  used  to  describe  these  processes  are  presented 
In  Appendix  I,  Section  A. 

A basic  characteristic  of  thermal  compressors  that  emerges  from  our 
analysis  Is  the  relationship  between  volumetric  efficiency  and  pressure 
ratio.  This  relationships  Is  shown  In  Figure  9.  For  a given  pressure 
ratio  and  flow  rate,  the  volumetric  efficiency  determines  cylinder  size. 
Hence,  this  relationship  Is  Important  In  sizing  the  device.  Figure  9 
shows  that  there  Is  a direct  trade-off  between  volumetric  efficiency  and 
pressure  ratio.  The  ambient  end  pressure-volume  diagrams  for  the  two 
extreme  cases  (one  In  which  the  volumetric  efficiency  Is  a maximum  and 
the  pressure  ratio  Is  unity,  and  the  other  In  which  the  pressure  ratio 
Is  maximum  and  the  volumetric  efficiency  Is  zero) , are  shown  In  Fig- 
ure 10.  Maximum  volumetric  efficiency  Is  achieved  when  the  gas  Is  not 
compressed,  l.e.,  when  the  ratio  of  discharge  pressure  to  suction 
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pressure  Is  one.  In  this  case  a large  portion  of  the  gas  entering  the 
cylinder  during  the  suction  stroke  Is  forced  out  of  the  cylinder  during 
the  discharge  stroke  at  slightly  above  suction  pressure.  Maximum  pres- 
sure ratio  Is  achieved  when  there  Is  no  through-flow,  l.e.,  when  the 
volumetric  efficiency  Is  zero.  In  this  case  the  gas  In  the  cylinder 
simply  shuttles  back  and  forth  between  the  hot  and  ambient  ends. 

Neither  of  these  extremes  are  useful  In  our  application.  Hence,  one  of 
the  major  design  choices  is  to  select  a pressure  ratio  per  cylinder 
that  yields  a reasonable  volumetric  efficiency  and  cylinder  size. 

In  Figure  9,  the  vertical  and  horizontal  Intercepts  are  determined 
from  Equation  A5  in  Appendix  I.  For  the  hot  end  temperature  of  1200#F 
and  an  ambient  temperature  of  80°F,  representative  values  of  these  two 
Intercepts  are  a vertical  Intercept  (l.e.,  the  maximum  volumetric 
efficiency)  of  60Z,  and  a horizontal  Intercept  (l.e.,  the  maximum  pres- 
sure ratio)  of  2.0.  A practical  design  point  will  lie  In  the  mid-range 
of  these  values,  l.e.,  a volumetric  efficiency  of  25  to  35%  and  a pres- 
sure ratio  of  1.4  to  1.6. 


We  have  considered  three  types  of  matrices  for  the  regenerator: 
packed  spheres,  stacked  screens,  and  cylindrical  tube  bundles.  The  heat 
transfer  and  pressure  drop  characteristics  of  these  matrices  are  repre- 
sented by  the  usual  correlations  of  friction  factor  and  Colburn  factor 
versus  Reynolds  number.  We  have  based  the  thermal  performance  of  the 
regenerator  on  the  mean  flow  rate  through  It  during  a cycle,  as  Is 
commonly  done  In  the  analysis  of  similar  devices.  However,  we  have 
calculated  the  pressure  drop  as  a function  of  displacer  position  to 
obtain  a more  complete  picture  of  how  the  pressure  force  on  the  dis- 
placer varies  throughout  the  cycle.  The  equations  describing  the  re- 
generator are  presented  In  Appendix  I,  Section  B. 


Heat  Exchangers 


We  have  considered  two  types  of  hot  and  ambient  heat  exchanger  con- 
figurations— annular  and  multipassage,  as  shown  in  Figure  11.  The  annular 
configuration  Is  simply  a thin  annular  passage,  while  the  multipassage 
configuration  consists  of  a number  of  small  tubes,  with  either  square 
or  circular  cross  section,  In  parallel.  The  thermal/f luld  processes  In 
the  heat  exchanger  are  analyzed  In  a manner  similar  to  the  regenerator, 
l.e.,  thermal  performance  Is  based  on  the  mean  flow  rate,  but  pressure 
drop  Is  calculated  as  a function  of  the  position  of  the  displacer. 

The  Reynolds  number  and  the  length-to-hydraulic  diameter  ratio  for 
these  exchangers  are  such  that  we  can  assume  that  fully  developed  lami- 
nar flow  will  exist  throughout  most  of  the  stroke. 


It  Is  worth  noting  that  the  heat  transferred  In  the  hot  and  ambient 
heat  exchangers  Is  much  less  than  the  heat  transferred  In  the  regenera- 
tor. Hence,  the  design  of  these  exchangers,  though  Important,  Is  not 
nearly  as  critical  as  that  of  the  regenerator.  The  relative  amount  of 
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FIGURE  1 1 HEAT  EXCHANGER  CONFIGURATIONS 


heat  transferred  in  the  heat  exchangers  and  in  the  regenerator  may  be 
seen  in  Figure  12,  which  is  a dlagramatic  representation  of  the  energy 
flows  in  the  compressor.  The  vertical  (temperature)  axis  on  this  figure 
is  approximately  to  scale,  so  vertical  distances  are  approximately  pro- 
portional to  heat  transferred.  From  this  figure,  one  can  see  that  the 
heat  transferred  in  the  regenerator  is  approximately  five  times  that 
transferred  in  the  ambient  heat  exchanger  and  the  hot  heat  exchanger. 

The  equations  used  in  analyzing  the  heat  exchangers  are  presented 
in  Appendix  I,  Section  C. 


5.  Displacer  Power  Input 

The  pressure  difference  between  the  hot  and  ambient  ends  acting  on 
the  area  A2,  constitutes  a force  thaf  is  always  opposite  in  direction  to 
the  velocity  of  the  displacer.  this  is  a damping  force  which 

removes  energy  from  the  displacer.  The  pressure  in  the  hot  end  acting  on 
the  projected  shaft  cross  section  area  (A3  - A2)  constitutes  a force 
that,  when  integrated  around  the  cycle,  provides  an  energy  input  to  the 
displacer.  The  difference  between  these  two  energy  values  determines  the 
net  power  input  which  is  required  to  sustain  motion  of  the  displacer.  A 
net  power  input  to  the  displacer  is  required  for  most  practical  designs. 
However,  in  extreme  cases,  in  which  either  the  pressure  drop  through  the 
regenerator  is  unusually  small,  or  the  area  (A3  - A2)  is  large,  it  is 
possible  to  obtain  a net  power  output  from  the  displacer. 

The  equations  used  to  determine  the  displacer  power  input  are  pre- 
sented in  Appendix  I,  Section  F. 

6.  Hot  End  Power  Input 

The  Ideal  heat  input  to  a thermal  compressor  is  the  same  as  the 
heat  supplied  to  a reversible  engine  which  is  driving  an  Isothermal  com- 
pressor. In  principle,  then,  the  thermal  compressor  is  a reversible 
device.  In  practice,  however,  a number  of  thermal  loss  mechanisms  pre- 
vent a close  approach  to  ideal  performance.  These  losses  are: 

a.  The  regenerator  loss  associated  with  enthalpy  flux  from 
the  hot  end  to  the  ambient  end. 

b.  Shuttle  heat  transfer 

c.  Displacer/cylinder  clearance  pumping  loss 

d.  Conduction  through  the  displacer 

e.  Conduction  through  cylinder  walls 

f.  Conduction  through  the  regenerator  shell (s) 

g.  Conduction  through  the  regenerator  matrix 
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h.  Conduction  through  insulation  surrounding  the  hot  end. 

The  equations  used  to  evaluate  these  elements  of  the  power  input 
are  presented  in  Appendix  I,  Section  E. 

C.  COMPUTER  PROGRAM 


Approximately  90  equations  are  required  to  define  the  thermo- 
dynamic processes  which  occur  in  a thermal  compressor,  and  to  design  the 
cylinder  assembly.  Further,  41  input  variables  must  be  specified  in 
order  to  perform  the  calculations.  In  order  to  expedite  calculations 
and  in  order  to  be  able  to  examine  the  effect  of  varying  a number  of  the 
more  important  variables,  the  governing  equations  were  programmed  for  a 
digital  computer.  The  program  performs  a straight  line  computation; 
starting  with  the  specified  input  parameters,  it  solves  all  the  equations 
in  sequence  and  prints  out  the  computational  results.  Since  there  is  no 
optimization  performed  by  the  program,  optimization  of  a machine  con- 
sists of  varying  input  parameters  over  a range,  examining  the  results, 
and  selecting  a design  which  satisfies  an  optimization  criterion  such  as 
minimum  thermal  power,  minimum  electrical  power,  etc.  The  computer  pro- 
gram is  described  in  detail  in  Appendix  II,  and  the  results  of  a number 
of  runs  performed  with  it  are  described  in  the  following  paragraphs. 


D.  RESULTS 


1.  Selection  of  the  Number  of  Cylinders 

The  computational  effort  in  the  parametric  analysis  was  directed 
specifically  toward  designing  a thermal  compressor  to  meet  a specified 
set  of  performance  requirements — those  determined  in  the  thermodynamic 
cycle  analysis.  The  first  step  in  this  process  was  to  select  the  number 
of  cylinders  for  the  compressor.  This  involved,  first,  selecting  the 
number  of  stages  of  compression  and  then  selecting  the  number  of  cylin- 
ders for  each  stage. 

In  determining  the  number  of  stages  of  compression,  one  must  con- 
sider the  flow/pressure  ratio  characteristic  for  a single  stage  of  a 
thermal  compressor.  Figure  9 showed  a plot  of  volumetric  efficiency 
versus  pressure  ratio.  On  this  plot,  a representative  value  for  the 
vertical  intercept  (i.e.,  the  maximum  volumetric  efficiency)  is  60%, 
and  a representative  value  of  the  horizontal  Intercept  (i.e.,  the  maxi- 
mum pressure  ratio)  is  2.0.  A practical  design  point  must,  of  course, 
lie  between  these  extremes.  Since  the  overall  compression  ratio  for  the 
compressor  is  2.0,  it  is  impractical  to  use  a single  stage  machine.  If 
two  stages  were  used,  each  stage  would  have  a compression  ratio  of  1.41 
(the  square  root  of  2);  and  if  three  stages  were  used,  each  stage  would 
have  a compresison  ratio  of  1.26  (the  cube  root  of  2).  The  correspon- 
ding volumetric  efficiencies  are  36%  and  44%.  Two  stages  of  compression 
were  selected  because  this  results  in  less  machinery  than  a machine  with 
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more  stages  and  because  the  displacer  dimensions  are  reasonable,  even 
though  the  volumetric  efficiency  Is  lower  than  for  a three  stage 
machine . 

The  next  question  which  must  be  resolved  Is  the  number  of  cylinders 
to  use  for  each  stage  and  the  overall  arrangement  of  these  cylinders  and 
their  drive  assemblies.  In  this  decision,  the  following  Items  must  be 
taken  Into  consideration: 

a.  The  desire  to  minimize  thermal  losses. 

b.  The  requirement  to  mount  displacers  In  pairs  on  a common 
center  line  to  achieve  dynamic  balance  of  the  machine. 

c.  The  desire  to  minimize  the  length  of  the  compressor 
assembly  for  ease  of  Integration  of  the  refrigerator  Into 
the  spacecraft. 

d.  The  desire  to  minimize  the  number  of  components  In  order 
to  maximize  the  reliability  of  the  refrigerator  and  mini- 
mize Its  size,  complexity,  and  cost. 

e.  The  desire  to  have  a small  displacer  diameter  In  order  to 
minimize  Its  fabrication  problems. 

Three  compressor  configurations  were  considered.  They  are: 

a.  Configuration  I,  which  utilizes  one  cylinder  per  stage 
(for  a total  of  two  cylinders)  with  a drive  assembly  for 
each  displacer.  This  results  In  two  compressor  modules, 
each  Incorporating  tme  cylinder  and  one  drive. 

b.  Configuration  II,  which  utilizes  two  cylinders  per  stage 
(for  a total  of  four  cylinders)  with  a double-ended  drive 
assembly  for  each  pair  of  displacers.  This  results  In 
two  compressor  modules,  each  Incorporating  two  cylinders 
and  one  drive. 

c.  Configuration  III,  which  utilizes  two  cylinders  per  stage 
(for  a total  of  four  cylinders)  with  a drive  assembly  for 
each  displacer.  This  results  In  four  compressor  modules, 
each  Incorporating  one  cylinder  and  one  drive. 

These  three  configurations  are  shown  schematically  In  Figure  13. 

The  three  configurations  are  compared  In  Table  VI.  Configuration  I 
contains  the  minimum  number  of  components,  and  results  In  the  shortest 
package  length.  It  also  permits  one  heat  source,  which  will  simplify 
the  hot  end  design  If  heat  Is  transferred  to  the  hot  cylinders  by,  for 
Instance,  a heat  pipe.  Additionally,  since  Configuration  I utilizes 
the  minimum  number  of  cylinders,  the  heat  losses  from  the  hot  end  will 
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TABLE  VI 


be  minimized  because  the  parasitic  heat  losses  from  a single  cylinder 
compressing  the  full  flow  of  gas  are  less  than  the  parasitic  heat  losses 
from  two  cylinders  each  of  which  compresses  half  of  the  full  flow.  The 
displacers  in  Configuration  I are  larger  than  the  displacers  in  Configu- 
rations II  and  III,  but  they  are  not  excessively  large  considering  pre- 
vious experience  (see  Reference  9).  Therefore,  since  Configuration  I is 
reasonable  in  size,  and  since  it  results  in  the  simplest,  lowest  loss 
design,  this  configuration  was  selected  for  the  compressor.  As  a result 
of  this  selection,  the  parametric  analysis  reduced  to  determining  the 
best  configuration  for  compressor  cylinders  to  meet  the  following  spec- 


ifications: 

First  Stage  Second  Stage 

Cylinder  Cylinder 


Mass  Flow  Rate,  lb /hr  3.66  3.66 

Inlet  Pressure,  psia  29.4  41.6 

Discharge  Pressure,  psia  41.6  58.8 


In  accordance  with  previous  experience  with  Vuilleumier  cycle 
refrigerators,  we  have  used  a hot  end  gas  temperature  of  1,200#F.  For 
purposes  of  design,  an  80°F  ambient  end  gas  temperature  was  used. 

2.  Preliminary  Screening  Runs 

Initially,  a number  of  screening  runs  on  the  first  stage  were  per- 
formed with  the  thermal  compressor  computer  program,  with  the  object  of 
determining  which  variables  are  most  critical,  which  variables  could  be 
set  early  on,  and  which  variables  would  have  to  be  examined  further.  A 
further  object  was  to  generate  enough  information  so  that  some  of  the 
input  variables  could  be  fixed.  In  these  runs  the  following  items  were 
varied:  running  speed,  the  ratio  of  stroke  to  piston  diameter,  regen- 

erator length,  and  hot  cylinder  length.  The  regenerator  length  was  the 
same  as  the  hot  cylinder  length  in  these  runs.  Initial  estimates  (based 
on  hand  calculations  and  experience)  were  used  for  the  other  input  var- 
iables, and  remained  fixed  during  these  runs.  The  insulation  heat  leak 
(normally  calculated  from  the  cylinder  geometry)  was  fixed  at  30  watts, 
a value  which  can  reasonably  be  achieved  in  practice. 

The  initial  runs  showed  that  there  are  two  major  sources  of  heat 
loss  from  the  hot  end,  shuttle  heat  transfer  and  regenerator  ineffec- 
tiveness. The  way  in  which  these  losses  vary  is  shown  in  Figure  14. 

The  trends  which  are  apparent  from  this  figure  are: 

a.  The  shuttle  loss  decreases  rapidly  with  increasing  speed 
and  less  rapidly  with  decreasing  stroke-to-d lame ter  ratio. 

b.  The  shuttle  loss  decreases  with  increasing  hot  cylinder 
length,  and  is  unacceptably  high  for  a 2 inch  hot  cylinder 
length  at  the  lower  end  of  the  speed  range  considered 
(300  CPM) . 
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FIGURE  14 


c.  The  regenerator  loss  Increases  gradually  with  Increasing 
speed,  and  is  unaffected  by  the  stroke-to-diameter  ratio. 


d.  The  regenerator  loss  decreases  with  increasing  regenerator 
length . 

Based  on  the  foregoing  results  from  the  first  stage,  we  concluded: 
1)  that  the  least  power  design  would  result  if  the  regenerator  length 
could  be  varied  independent  of  the  hot  cylinder  length,  and  2)  that  the 
hot  cylinder  length  should  be  as  long  as  possible.  As  a result,  the 
hot  cylinder  length  was  set  at  6 inches  in  all  subsequent  analysis. 

After  fixing  the  hot  cylinder  length,  we  made  further  runs  in  which 
the  speed  and  the  stroke-to-diameter  ratio  were  varied. 

Figure  15  shows  the  results  of  one  set  of  these  runs.  It  is  a plot 
of  displacer  shaft  power,  displacer  diameter,  and  total  hot  end  power  as 
a function  of  speed  for  the  case  where  the  regenerator  length  is 
2 inches,  the  hot  cylinder  length  is  6 inches,  and  the  intermediate  void 
fraction  is  0.5.  Data  are  plotted  for  two  values  of  stroke-to-diameter 
ratio.  Similar  plots  for  other  values  of  regenerator  length  and  inter- 
mediate void  fraction  show  the  same  trends,  so  Figure  15  may  be  con- 
sidered as  representative.  Several  trends  are  apparent  from  this 
figure. 

a.  The  shaft  power  increases  with  increasing  speed,  while  the 
hot  end  power  decreases  with  increasing  speed.  Thus,  in 
selecting  a speed,  a balance  must  be  struck  between  these 
two  factors. 

b.  Displacer  diameter 'decreases  with  increasing  speed. 

c.  Thermal  power  is  not  a strong  function  of  the  stroke-to- 
diameter  ratio,  but  the  displacer  diameter  is  smaller  at 
higher  values  of  the  stroke-to-diameter  ratio. 

From  these  results,  we  concluded  that  the  most  compact  design  would 
result  if  the  stroke-to-diameter  ratio  were  at  the  upper  end  of  the 
range  considered,  so  this  ratio  was  set  at  0.2.  Additionally,  we  con- 
cluded that  running  speeds  above  900  CPM  would  result  in  unacceptably 
high  shaft  power  requirements,  so  we  narrowed  the  range  of  speeds  to  the 
range  300  CPM-900  CPM.  Having  set  these  parameters,  we  then  proceeded 
to  examine  the  effect  of  regenerator  design  parameters. 

3.  Regenerator  Analysis 

Since  the  regenerator  has  a major  influence  on  the  displacer  shaft 
power  and  the  hot  end  power,  we  made  a number  of  computer  runs  in  order 
that  we  might  select  a regenerator  design  which  would  result  in  a 
reasonable  balance  between  these  two  input  powers.  In  these  runs,  the 
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parameters  which  determine  regenerator  performance  were  varied.  These 
parameters  are:  speed,  intermediate  void  volume,  regenerator  matrix 

size,  and  regenerator  length.  All  runs  were  made  with  a packed  sphere 
matrix,  since  this  matrix  appears  to  be  preferable  to  screens.  (The 
pressure  drop  and  heat  transfer  characteristics  of  packed  spheres  are 
essentially  the  same  as  for  screens,  but  the  void  fraction  is  substan- 
tially less — .37  to  .39  for  spheres  compared  to  .60  for  screens.) 

Results  of  a representative  series  of  runs  are  shown  in  Figure  16,  in 
which  hot  end  power  and  displacer  shaft  power  are  plotted  against  sphere 
diameter  for  different  values  of  regenerator  length.  The  data  shown  are 
for  a speed  of  600  CPM  and  an  intermediate  void  fraction  of  0.5.  It  can 
be  seen  that  for  a given  regenerator  length,  as  the  sphere  diameter 
increases,  the  shaft  power  decreases  while  the  hot  end  power  increases. 
Also,  lengthening  the  regenerator  at  a fixed  matrix  size  increases  the 
displacer  shaft  power  and  decreases  the  hot  end  power.  These  data 
reconfirm  a basic  characteristic  of  the  thermal  compressor — there  is  a 
tradeoff  between  displacer  shaft  power  and  hot  end  power,  with  a reduc- 
tion in  one  form  of  power  being  accompanied  by  an  increase  in  the  other. 


For  each  value  of  displacer  shaft  power,  there  is  a regenerator 
packing  diameter  and  a regenerator  length  which  result  in  a minimum 
hot  end  power.  The  dotted  lines  on  the  lower  plot  of  Figure  16  show 
this  minimum.  These  dotted  lines  are  simply  lines  of  constant  displacer 
shaft  power  transferred  from  the  upper  plot  on  the  figure.  Note  that 
the  optimum  regenerator  matrix  (i.e.,  the  one  which  results  in  minimum 
hot  end  power)  depends  upon  the  value  of  displacer  shaft  power  chosen. 

If  one  desires  a low  displacer  shaft  power,  then  one  must  employ  a 
relatively  short  regenerator  packed  with  relatively  large  spheres,  and 
must  accept  a relatively  high  hot  end  power.  On  the  other  hand,  if  one 
is  willing  to  accept  a higher  displacer  shaft  power,  one  can  employ  a 
relatively  long  regenerator  packed  with  relatively  small  spheres,  and 
will  achieve  a relatively  low  hot  end  power.  Plots  similar  to  Figure  16 
were  made  for  speeds  of  300,  600,  and  900  CPM  at  two  values  of  inter- 
mediate void  fraction,  0.3  and  0.5.  From  these  plots,  the  minimum  hot 
end  power  was  determined  as  a function  of  displacer  drive  power.  The 
results  of  these  determinations  are  shown  on  Figure  17,  on  which  hot  end 
power  is  plotted  against  displacer  shaft  power  at  several  values  of 
speed.  Two  sets  of  curves  are  shown- — one  for  each  of  two  values  of 
intermediate  void  fraction.  Several  things  are  apparent  from  Figure  17: 


a. 


b. 


Hot  end  power  decreases  as  displacer  power  increases. 

At  all  values  of  displacer  shaft  power,  there  is  a speed 
which  results  in  minimum  hot  end  power.  This  speed  is 
600  CPM  over  most  of  the  range  of  displacer  drive  power. 


At  a speed  of  600  CPM,  the  hot  end  power  is  virtually 
independent  of  intermediate  void  fraction. 
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The  Baseline  Design  Point  was  selected  primarily  by  considering  the 
information  contained  in  Figure  17 , and  the  conclusions  resulting  from 
it.  A speed  of  600  CPM  was  selected,  since  this  results  in  minimum  hot 
end  power.  At  a speed  of  600  CPM,  an  intermediate  void  fraction  of  0.5 
results  in  somewhat  lower  hot  end  power,  so  an  intermediate  void  fraction 
of  0.5  was  selected. 


The  selection  of  regenerator  matrix  si?e  and  regenerator  length  is 
not  so  straightforward,  since  these  parameters  determine  the  displacer 
shaft  power  and  the  hot  end  power  which  will  be  required  to  run  the 
compressor.  Thus,  implicit  in  the  selection  of  the  regenerator  length 
and  matrix  size  is  a value  judgment  regarding  the  relative  importance  of 
hot  end  power  (thermal  power)  and  displacer  shaft  power  (electrical 
power).  A prime  feature  of  the  thermal  compressor  is  that  it  reduces 
the  electric:^  drive  power  requiiad  for  the  refrigerator  by  substituting 
thermal  powAr^or  some  of  the  electric  power  required  by  an  electric 
compressor.'  As  noted  previously  low  values  of  displacer  drive  power 
are  accompanied  by  high  values  of  hot  end  power  and  vice  versa.  Thus, 
in  selecting  a baseline  design  point,  we  must  strike  a balance  between 
the  electric  power  and  the  thermal  power.  In  seeking  this  balance,  we 
ask  the  following  two  questions:  1)  is  there  a displacer  drive  power 

below  which  thermal  power  increases  rapidly  and  above  which  thermal 
power  decreases  rather  slowly?;  and  2)  considering  the  refrigeration 
system  as  a whole,  is  the  electric  input  power  required  to  overcome 
pressure  differences  across  the  two  displacers  in  proportion  to  the 
electric  power  required  for  other  portions  of  the  refrigerator  system? 

A nominal  value  of  15  watts  of  drive  power  per  displacer  (or  30  watts 
for  the  two  displacers)  seems  tc  satisfy  both  these  criteria  reasonably 
well.  As  to  the  first  question.,  on  Figure  17  one  can  see  that,  for  a 
speed  of  600  CPM  and  an  intermediate  void  fraction  of  0.5,  the  hot  end 
power  increases  rapidly  at  displacer  drive  powers  of  less  than  approxi- 
mately 15  watts,  but  that  it  increases  rather  slowly  at  displacer  drive 
powers  greater  than  approximately  15  watts.  As  regards  overall  electric 
power,  the  entire  thermal  compressor  refrigerator  is  estimated  to  re- 
quire 92  watts  of  electric  power  (see  Section  VI),  20  of  which  are 
required  by  the  expander  assembly,  and  thus  are  completely  independent 
of  the  compressor  design  It  does  not  seem  out  of  proportion  to  allot 
30  watts  (approximately  one  third  of  the  total  power)  to  overcome  pres- 
sure drop  losses  in  the  compressor  assembly.  Based  on  the  foregoing, 
we  adopted  a target  of  15  watts  for  a displacer  drive  power  for  each 
compressor  stage.  During  the  compressor  design,  we  selected  a regenera- 
tor length  and  matrix  size  which  corresponded  to  this  design  point. 

For  the  first  stage,  this  results  in  a regenerator  length  of  two 
inches  and  a matrix  comprised  of  .040  inch  diameter  spheres  (from 
Figure  16). 
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Summarizing,  then,  a number  of  variables  were  selected  and  fixed 
during  the  course  of  the  parametric  analysis.  The  fixing  of  these 
variables  specifies  a compressor  we  have  termed  the  baseline  design. 

The  primary  variables  for  this  design  (that  is,  the  ones  which  have  been 
discussed  in  this  section)  are  repeated  below. 


Number  of  Stages 

2 

Number  of  Cylinders  Per  Stage 

1 

Stroke-to-Dlameter  Ratio 

.2 

Hot  Cylinder  Length 

6 In 

Hot  End  Gas  Temperature 

1200*F 

Ambient  Gas  Temperature 

80#F 

Speed 

600  CPM 

Intermediate  Void  Fraction 

0.5 

Displacer  Shaft  Power 

15  watts 

Other  variables  which  were  fixed  during  the  course  of  the  detailed  com- 
pressor design  are  discussed  in  the  next  section  and/or  are  displayed  in 
the  sample  problems  in  Appendix  II. 


The  baseline  design  compressor  was  examined  in  somewhat  more  detail 
in  the  course  of  actually  laying  out  the  design  of  the  compressor.  The 
examination  took  the  form  of  making  final  adjustments  to  hot  and  ambient 
void  fractions,  and  in  taking  into  account  in  more  detail  the  pressure 
drop  in  gas  passages  which  interconnect  the  working  volumes,  the  regen- 
erators, and  the  two  heat  exchangers.  As  a result  of  these  further 
refinements  of  the  design,  the  actual  compressor  design  performance  is 
somewhat  different  than  shown  in  the  parametric  curves  in  this  section. 
The  final  baseline  design  is  described  in  the  following  section. 
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V.  DESIGN  OF  BASELINE  THERMAL  COMPRESSOR 


A.  GENERAL 


We  have  generated  the  preliminary  design  of  a rotary-reciprocating 
thermal  compressor  to  meet  performance  specifications  listed  in  Table 
II;  that  is,  a compressor  to  compress  3.66  lb/hr  of  helium  gas  at  300  K 
from  29.4  lb /in2  inlet  pressure  to  58.8  lb /in2  discharge  pressure.  The 
design  of  the  compression  cylinders  has  utilized  the  results  of  the 
parametric  analysis  described  in  the  preceding  section.  The  major 
emphasis  in  the  design  of  the  cylinder  assemblies  has  been  on  the  dis- 
placers, heat  exchangers,  regenerators,  and  valving.  Relatively  little 
emphasis  has  been  directed  toward  the  type  of  heat  source  used  or  its 
integration  with  the  cylinders.  The  design  of  the  drive  assemblies  has 
drawn  upon  our  previous  experience  with  rotary— reciprocating  compressors 
and  expanders.  The  drive  assemblies  have  been  designed  along  the  lines 
of  previous  devices — that  is,  they  utilize  a rotary  motor  and  linear 
actuator  to  drive  the  displacer,  and  contain  gas  springs  to  provide  a 
restoring  force  to  the  displacer. 

The  first  stage  of  the  baseline  design  thermal  compressor  is  shown 
in  Figure  18,  and  the  second  stage  is  shown  in  Figure  19.  The  compres- 
sor is  a two-cylinder  machine,  with  the  two  cylinders  arranged  in  line 
on  a common  axis.  The  hot  ends  of  both  cylinders  are  in  the  center  of 
the  assembly,  and  are  surrounded  by  a vacuum  shell  containing  high  tem- 
perature insulation.  Electrical  cables  connect  the  drive  assemblies  to 
a power  conditioning  equipment  (PCE)  assembly.  The  PCE  assembly  con- 
trols the  power  to  the  linear  actuators  and  rotary  motors  in  the  drive 
assembly.  The  characteristics  of  the  baseline  design  are  summarized  in 
Table  VII.  The  following  paragraphs  describe  the  various  elements  in 
the  compressor. 


FIRST  STAGE  CYLINDER  ASSEMBLY 


The  first  stage  cylinder  assembly  consists  of  the  cylinder,  the 
displacer,  the  regenerator,  the  hot  and  ambient  heat  exchangers,  and  a 
porting  system.  The  pressure-containing  shell  of  the  cylinder  bolts 
onto  a large  flange  on  the  end  of  the  vacuum  shell.  There  is  a .010 
inch  thick  cylindrical  section  in  this  piece  immediately  adjacent  to 
the  attachment  flange.  This  section  is  the  outer  shell  of  the  regener- 
ator. The  wall  is  thin  in  this  zone  to  minimize  axial  conduction  down 
the  wall  of  the  regenerator,  since  the  entire  temperature  gradient  in 
the  machine  appears  across  this  section.  The  wall  thickness  of  the 
cylinder  is  Increased  in  the  zone  between  the  regenerator  and  the  head 
of  the  cylinder.  This  zone  contains  two  gas  flow  passages  and  the 
hot  heat  exchanger.  The  head  of  the  cylinder  is  also  thick—both  to 
reduce  stresses  and  to  provide  a good  conductive  heat  transfer  path 
across  the  entire  cylinder  head.  A hollow  tubular  insert  in  the  shell 
of  the  cylinder  contains  the  cylinder  itself  and  the  inner  boundaries 
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TABLE  VII 


SUMMARY  OF  BASELINE  DESIGN  THERMAL  COMPRESSOR 


Process  Fluid 


Flow  Rate 


Hot  End  Gas  Temperature 


Ambient  End  Gas  Temperature 


First  Stage 


Inlet  Temperature  °F 


Inlet  Pressure,  LB/IN2 


Discharge  Pressure,  LB/IN 


Displacer  Diameter,  IN 


Stroke,  IN 


Reciprocating  Speed,  Cycles /Min 


Rotational  Speed,  Rev/Min 


Input  Power  (Total  for  Both  Stages) 


Thermal,  Watts 


Electrical,  Watts 


Overall  Diameter v * , IN 


Overall  Length v , IN 


Weight  Pounds 


'Not  Including  power  conditioning  equipment. 


1,064 


r ' ' ' , 


Helium  Gas 


3.66  LB/HR 


1,200*F 


Second  Stage 


P9PhP0S 


of  the  regenerator,  gas  flow  passages  and  the  hot  heat  exchanger.  This 
Insert  contains  a 6 Inch  long,  .010  Inch  thick  section  adjacent  to  the 
displacer,  over  which  the  entire  hot  end  temperature  gradient  appears. 

It  contains  a shorter,  .010  Inch  thick  section  In  the  vicinity  of  the 
regenerator.  The  hollow  annulus  In  the  Insert  Is  evacuated  and  filled 
with  insulation  to  reduce  heat  transfer  from  the  cylinder  wall  to  the 
hot  heat  exchanger  and  regenerator.  This  is  necessary  since  the 
temperature  gradient  in  the  cylinder  wall  is  different  than  the  temper- 
ature gradient  in  the  hot  heat  exchanger  and  the  regenerator.  The 
outer  cylinder  and  the  Insert  are  both  made  of  Inconel  718. 

The  displacer  consists  of  two  pieces,  an  aluminum  disc  at  the  am- 
bient end  and  an  Inconel  718  hat  section  which  is  mounted  to  the  alumi- 
num disc.  The  displacer  Is  5.6  inches  in  diameter.  It  Is  hollow,  and 
is  evacuated  and  filled  with  Insulation.  The  wall  of  the  displacer  Is 
reduced  to  the  extent  possible  considering  that  It  Is  a vacuum  vessel. 

It  is  .030  inches  thick  and  contains  several  Internal  stiffening  rings. 
There  is  a .020  inch  radial  gap  between  the  displacer  and  the  cylinder 
in  the  zone  over  which  the  temperature  gradient  occurs  to  reduce  the 
shuttle  heat  transfer  losses.  This  gap  is  also  sufficient  to  eliminate 
any  mechanical  problem  associated  with  distortion  of  the  hot  ends  of 
the  displacer  or  the  cylinder.  Bypass  flow  through  the  annular  gap 
between  the  displacer  and  the  cylinder  is  limited  by  a clearance  seal 
at  the  ambient  end  of  the  displacer.  The  radial  gap  in  this  seal  area 
is  .001  inches.  The  displacer  is  aligned  in  the  cylinder  by  the  two  gas 
bearings  in  the  drive  assembly.  The  clearance  in  these  bearings  is  less 
than  the  clearance  in  the  displacer  seal,  and  concentricity  tolerances 
are  such  that  there  is  no  contact  between  displacer  and  the  cylinder. 

The  regenerator  is  annular  in  shape.  It  is  two  inches  long,  and 
has  an  inside  diameter  of  6.6  inches  and  an  outside  diameter  of  7.9 
inches.  It  is  packed  with  0.40  inch  diameter  stainless  steel  spheres 
which  are  spring  loaded  in  the  axial  direction  to  prevent  the  packing 
from  rattling.  The  effectiveness  of  the  regenerator  is  calculated  to 
be  97%. 

The  hot  heat  exchanger  is  an  annulus  defined  by  the  inside  of  the 
outer  cylinder  wall  and  the  outside  of  the  tubular  insert.  It  is 
1.3  inches  long  and  has  an  annular  gap  of  .055  inches.  Its  effective- 
ness is  specified  to  be  90%. 

The  ambient  heat  exchanger  also  has  an  annular  configuration.  It 
is  fabricated  by  joining  (brazing  or  welding)  a tapered  insert  into 
the  housing  prior  to  final  machining  of  the  housing.  The  gap  between 
the  insert  and  the  main  housing  serves  as  the  ambient  heat  exchanger. 
This  heat  exchanger  is  1.8  inches  long  and  has  a radial  gap  of  .035 
inches.  Its  effectiveness  is  specified  to  be  90%. 
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The  porting  Is  contained  In  the  displacer  drive  shaft . It  consists 
of  two  transverse  holes  which  are  connected  by  an  axial  hole  on  the 
shaft  center  line.  The  transverse  hole  adjacent  to  the  displacer  Is 
always  In  communication  with  the  ambient  working  volume,  while  the  other 
transverse  hole  communicates  with  the  Inlet  and  discharge  ports  in  the 
housing  at  appropriate  axial  positions  for  correct  valving  action. 

Figure  20  Indicates  the  port  opening  and  port  closing  points.  It  shows 
the  pressure  in  the  first  and  second  stage  ambient  working  volumes  as  a 
function  of  displacer  position  for  the  condition  of  no  pressure  drop  In 
the  porting.  For  this  condition,  the  first  stage  inlet  ports  start  to 
open  .51  Inches  from  top  dead  center  and  are  completely  closed  at  bottom 
dead  center,  while  the  first  stage  discharge  ports  start  to  open  at  .43 
Inches  from  top  dead  center  and  are  completely  closed  at  top  dead  cen- 
ter. The  actual  location  of  the  two  port  opening  points  must  take  Into 
account  the  pressure  drop  through  the  ports  as  they  start  to  open.  When 
this  pressure  drop  Is  taken  into  account,  the  opening  point  of  the  suc- 
tion port  will  be  shifted  slightly  toward  top  dead  center,  and  the 
opening  point  of  the  discharge  port  will  be  shifted  slightly  toward 
bottom  dead  center.  There  will  be  little  or  no  change  In  the  location 
of  the  port  closing  points,  however,  since  the  piston  velocity,  and 
hence  porting  pressure  drop.  Is  very  low  at  the  two  extremes  of  stroke. 

The  porting  zone  Is  bounded  by  two  circumferential  grooves  In  the 
housing;  the  one  adjacent  to  the  Inlet  port  Is  communicated  to  Inlet 
pressure,  and  the  one  adjacent  to  the  discharge  port  Is  communicated 
to  discharge  pressure.  These  grooves  confine  the  axial  and  circumfer- 
ential pressure  gradients  in  the  porting  zone,  and  provide  a constant 
pressure  boundary  for  the  clearance  seal  on  one  side  and  the  gas 
bearing  on  the  other. 

The  thermal  input  power  to  the  hot  end  of  the  cylinder  must  be 
sufficient  to  provide  the  Ideal  P-V  energy  to  the  hot  working  volume  and 
also  to  make  up  for  thermal  losses.  For  the  first  stage,  the  ideal 
Input  power  Is  172  watts  and  the  losses  are  380  watts  (approximately 
2.2  times  the  Ideal  Input  power)  so  the  total  heat  input  to  the  hot  end 
is  552  watts.  A breakdown  of  this  power  Is  shown  In  Table  VIII.  Note 
that  the  major  losses  are  due  to  shuttle  heat  transfer  and  regenerator 
Ineffectiveness. 

In  Figure  18  the  heat  source  at  the  hot  end  of  the  cylinder  has 
been  shown  as  an  electric  heater.  This  form  of  heat  Input  has  been 
shown  for  Illustrative  purposes  only.  Other  heat  sources  could  also  be 
used  for  Instance,  focused  sunlight  or  heat  from  a radioisotope.  The 
type  of  heat  source  which  Is  most  appropriate  to  a specific  application 
and  the  means  of  transferring  this  heat  Into  the  hot  cylinder  (for 
Instance,  via  a heat  pipe  or  via  radiation)  would  depend  upon  system 
considerations.  Of  course,  the  thermal  Interface  at  the  hot  end  would 
have  to  be  tailored  to  the  heat  source  and  the  heat  transfer  means. 

The  fact  that  both  hot  cylinders  are  adjacent  to  each  other  will  sim- 
plify Integration  of  the  heat  source. 
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TABLE  VIII 

BREAKDOWN  OF  INPUT  POWER 
FOR  BASELINE  DESIGN  THERMAL  COMPRESSOR 
(Watts) 

Thermal  Poorer 


First  Stage 

Second  Stage 

Ideal  Input  Power 

172 

172 

* 

Losses: 

Shuttle  Heat  Transfer 

141 

118 

Displacer  Conduction 

22 

19 

Cylinder  Wall  Conduction 

7 

6 

Insulation  Conduction 

30 

30 

Pumping  Loss 

1 

2 

Regenerator  Heat  Loss 

111 

109 

Regenerator  Shell  Conduc- 
tion 

57 

49 

Regenerator  Matrix  Conduc- 
tion 

11 

7 

TOTAL  Losses 

380 

340 

TOTAL  Thermal  Power 

552 

512 

Displacer  Drive  Power 

First  Stage 

Second  Stage 

Pressure  Drop  Loss 

26 

25 

Input  from  Hot  End 

11 

10 

Net  Drive  Power 

15 

15 

I 


hkk 
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Table  VIII  also  gives  a breakdown  of  the  displacer  drive  power, 
and  indicates  that  the  26  watt  pressure  drop  loss  is  partially  supplied 
by  11  watts  of  power  from  the  hot  working  volume  (due  to  the  difference 
in  area  between  the  hot  working  volume  and  the  cold  working  volume). 

Thus,  15  watts  of  power  must  be  supplied  to  the  displacer  from  the 
drive  assembly. 

Heat  rejection  at  the  ambient  end  of  the  cylinder  is  accomplished 
by  circulating  a coolant  through  machined  passages  in  the  cylinder 
housing.  An  estimate  of  the  power  consumption  of  the  pump  in  the  heat 
rejection  system  has  been  included  in  Section  VII,  in  which  a refrigeration 
system  incorporating  a thermal  compressor  and  a refrigeration  system 
incorporating  an  electrical  compressor  are  compared. 


C.  SECOND  STAGE  CYLINDER  ASSEMBLY 

In  its  overall  configuration,  the  second  stage  cylinder  assembly 
(Figure  19)  is  the  same  as  the  first  stage  assembly;  the  only  differences 
being  that  some  of  the  dimensions  are  different . The  stroke  and  speed 
of  the  second  stage  are  the  same  as  the  first  «tage  (1.12  inches  and 
600  CPM,  respectively)  to  permit  dynamic  balance  of  the  compressor  as 
a whole.  A displacer  diameter  of  4.7  inches  is  required  to  give  the 
proper  throughput  at  this  speed  and  stroke.  The  annular  regenerator  is 
2 inches  long  and  has  an  inside  diameter  of  5.7  Inches  and  an  outside 
diameter  of  6.8  inches.  It  is  packed  with  .034  inch  diameter  stainless 
steel  spheres  and  its  effectiveness  is  calculated  to  be  97%.  The  hot  and 
ambient  heat  exchangers  are  of  an  annular  configuration.  The  hot  heat 
exchanger  is  1.0  inches  long  with  a radial  gap  of  .038  Inches  and  has  an 
effectiveness  of  90%,  while  the  ambient  heat  exchanger  is  1.7  inches 
long,  with  a radial  gap  of  .029  inches  and  has  an  effectiveness  of  90%. 

A pressure  versus  displacement  diagram  and  the  porting  dimensions  were 
shown  in  Figure  20.  It  can  be  seen  from  this  figure,  for  the  condition 
of  no  pressure  drop  in  the  porting,  the  inlet  port  starts  to  open 
.51  Inches  from  top  dead  center  and  the  discharge  port  starts  to  open 
.43  inches  from  top  dead  center.  The  thermal  input  power  to  the  second 
stage  is  512  watts.  The  breakdown  of  this  power  was  shown  in  Table 
VIII.  As  with  the  first  stage,  the  major  losses  are  due  to  shuttle  heat 
transfer  and  regenerator  ineffectiveness.  The  displacer  drive  power  for 
the  second  stage  is  15  watts.  A breakdown  of  this  power  was  also  shown 
in  Table  VIII. 


D.  DRIVE  ASSEMBLIES 

The  two  drive  assemblies  are  identical.  They  each  contain  two  gas 
bearings  to  support  the  displacer,  a gas  spring  to  resonate  the  displacer 
mass,  a linear  actuator  to  drive  the  displacer,  a rotary  motor  to  rotate 
the  displacer,  and  an  LVDT  to  provide  an  axial  position  signal  to  the 
amplitude  control  circuitry.  The  material  of  construction  is  aluminum 
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throughout,  with  the  exception  that  the  linear  actuator  and  rotary 
motor  use  Vanadium  Permendur  in  the  magnetic  circuits. 

The  gas  bearings  are  360°  journal  bearings  which  are  one  inch  in 
diameter  by  two  Inches  long.  They  operate  at  the  first  stage  discharge 
pressure.  Each  bearing  has  a radial  clearance  of  75  mlcrolnches  and 
has  a load  carrying  capacity  of  12  pounds  at  an  operating  eccentricity 
ratio  of  0.3.  This  load  capacity  Is  sufficient  to  support  weight  loads, 
the  pneumatic  load  In  the  porting  area,  plus  magnetic  forces  from  the 
linear  actuator  and  rotary  motor.  All  surfaces  on  the  bearing  diameters 
(that  Is,  bearings,  seals,  and  porting)  are  coated  with  plasma-sprayed 
chromium  oxide  to  provide  wear  resistance  at  startup  and  shutdown. 

The  clearance  seals  between  the  bearings  and  the  gas  spring  cavi- 
ties limit  the  leakage  flows  from  the  gas  spring  cavities.  The  clear- 
ance seal  between  the  outboard  bearing  and  the  motor  cavity  limits  the 
leakage  from  the  bearing  (which  operates  at  first  stage  discharge  pres- 
sure) and  the  motor  cavity  (which  operates  at  the  Inlet  pressure  to  the 
stage).  All  seals  have  the  same  clearance  as  their  adjacent  bearings. 

The  gas  spring  must  have  an  equivalent  spring  constant  such  that 
the  resonant  frequency  of  the  6.8  pounds  displacer  Is  600  CPM — that  Is, 
an  equivalent  spring  constant  of  69  pounds  per  Inch.  For  the  case  In 
which  the  gas  spring  pressure  Is  referenced  to  the  first  stage  Inlet 
pressure  at  the  neutral  position,  this  equivalent  spring  constant  will 
result  from  a gas  spring  which  has  an  outside  diameter  of  1.7  Inches 
and  an  Inside  diameter  of  1.0  Inches  (the  bearing  diameter),  and  a dead 
volume  at  neutral  position  which  Is  four  times  one  half  of  the  displace- 
ment volume  of  the  gas  spring.  For  this  design,  the  pressure  In  the  gas 
spring,  cavity  will  fluctuate  between  20.2  psla  and  47.3  psla.  The 
power  loss  In  the  gas  spring  due  to  heat  transfer  effects  is  estimated 
to  be  In  the  range  of  5-10  watts  for  this  design.  To  be  conservative, 
the  larger  figure  has  been  used  in  determining  power  requirements.  Each 
gas  spring  cavity  Is  referenced  to  the  first  stage  inlet  pressure  each 
time  the  displacer  passes  through  the  midstroke  position.  This  is 
accomplished  by  a porting  system  located  on  the  large  diameter  of  the 
gas  spring.  There  are  also  clearance  seals  on  this  diameter  to  limit 
leakage  flow  between  the  gas  spring  porting  system  and  the  spring  cavi- 
ties. 

The  linear  actuator  provides  the  driving  power  to  reciprocate  the 
displacer.  It  drives  the  displacer  toward  the  hot  end  during  half  of 
the  stroke,  and  gas  pressure  forces  drive  the  displacer  toward  the 
ambient  end  during  the  other  half  of  the  stroke.  The  power  which  the 
linear  actuator  must  supply  to  the  displacer  Is  the  sum  of  all  the 
losses  minus  the  power  delivered  by  the  hot  working  volume.  This  power 
Is  calculated  to  be  25  watts  for  each  compressor  stage  (see  Table  IX) . 
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TABLE  IX 


ACTUATOR  DRIVE  POWER  REQUIREMENTS 


(Watts) 


First  Stage 


Second  Stage 


Pressure  Drop  across  Displacer 


Gas  Spring  Losses 


Viscous  Shear  In  Bearings  and 
Seals 


TOTAL  Losses 


Net  Actuator  Shaft  Power  Requirement 


Since  all  displacer  losses  are  viscous  in  nature,  the  peak  actuator 
force  will  be  required  at  the  point  in  the  reciprocating  cycle  where  the 
velocity  is  the  greatest — that  is,  at  midstroke.  This  peak  force  is 
14.1  pounds.  Iron  and  copper  losses  in  the  actuators  are  estimated  to 
be  5 watts,  so  the  electrical  input  power  to  each  actuator  is  30  watts 
(the  shaft  power  plus  the  actuator  losses).  The  actuator  design  is 
similar  to  actuators  which  have  been  built  in  the  past.  They  will  con- 
sist of  a large  number  of  tapered  laminations  bonded  into  the  ring. 

The  magnetic  material  in  the  actuators  is  Vanadium  Permendur. 

The  rotary  motor  must  have  sufficient  starting  torque  to  overcome 
static  friction  and  sufficient  running  torque  to  overcome  viscous  shear 
losses  in  the  bearings  and  seals.  It  must  be  able  to  develop  1.6  pound- 
inches  of  starting  torque  and  .018  pound-inches  of  running  torque.  At 


300  rpm,  the  rotational  speed  of  the  displacers,  the  running  torque  is 
equivalent  to  .01  watts  of  input  power.  A value  of  1 watt  of  input 
power  to  each  rotary  motor  has  been  used  in  determining  input  power 
requirements.  The  mechanical  design  of  the  rotary  motor  is  similar  to 
previous  devices.  It  will  incorporate  a pe’-f.v.nent  magnet  motor  and  a 
two-pole-pair  stator  of  conventional  construction.  The  magnetic  mate- 
rial will  be  Vanadium  Permendur. 


SUMMARY  OP  THE  DESIGN 


The  salient  features  of  the  design  of  the  baseline  thermal  com- 
pressor were  given  in  Table  VII.  The  breakdown  of  the  electrical  input 
power  and  the  weight  axe  given  in  Tables  X and  XI. 


TABLE  XI 

WEIGHT  BREAKDOWN  OF  BASELINE  THERMAL  COMPRESSOR 

(Pounds) 


Two  Moving  Masses 
Vacuum  Shell  and  Insulation 
Two  Cylinder  Assemblies 
Two  Housing  Assemblies^) 
Two  Linear  Actuator  Stators 
Two  LVDT's 

Two  Rotary  Motor  Stators 
Auxiliaries^) 


TOTAL 


^) Includes  Porting,  Gas  Spring,  and  Motor  Housings 


'Piping,  Valves,  Surge  Chambers,  Electrical  Fittings 


VI.  SYSTEM  COMPARISONS 


In  the  preceding  section,  we  discussed  the  baseline  thermal  compres- 
sor as  a system  component.  In  this  section,  we  compare  the  overall 
characteristics  of  a complete  refrigeration  system  which  Incorporates  a 
thermal  compressor  against  a system  which  Incorporates  an  electric  compres- 
sor. Both  systems  are  sized  to  meet  the  refrigeration  specifications  which 
were  listed  In  Table  I. 

In  Table  XII,  the  major  features  of  the  two  refrigeration  systems  are 
compared.  In  this  comparison,  the  refrigerator  Is  considered  to  consist 
of  the  compressor  and  expander  assemblies,  their  power  conditioning  equip- 
ment, and  the  Interconnecting  piping  between  the  two  assemblies.  It  does 
not  consider  either  the  weight  or  power  consumption  of  the  heat  rejection 
system,  or  the  weight  of  the  electrical  or  thermal  power  supplies. 

It  is  seen  that  the  system  with  a thermal  compressor  requires  consider- 
ably more  power  than  the  system  with  an  electric  compressor;  1,156  watts 
compared  to  587  watts,  or  twice  as  much  power  for  the  thermal  compressor 
system.  However  only  92  watts,  or  8%,  of  the  power  for  the  thermal  compres- 
sor system  is  in  the  form  of  electrical  energy.  This  is  15%  of  the  electri- 
cal power  required  by  the  electric  compressor  system.  From  these  power 
figures,  and  from  the  power  consumption  of  other  systems  compared  during 
the  course  of  this  study,  we  can  make  the  following  generalizations: 

1.  The  total  power  requirement  of  a thermal  compressor  refri- 
gerator Is  2-3  times  that  of  an  equivalent  electric 
compressor  refrigerator. 

2.  The  electrical  power  requirement  of  a thermal  compressor 
refrigerator  Is  15-25%  of  that  of  an  equivalent  electric 
compressor  refrigerator. 

3.  The  thermal  power  requirement  of  a thermal  compressor 
refrigerator  Is  2-2  1/2  times  the  electrical  power 
requirement  of  an  equivalent  electric  compressor 
refrigerator. 

The  size  of  the  two  refrigerators  are  virtually  the  same,  while  the 
weight  of  the  thermal  compressor  refrigerator  Is  somewhat  greater  (approxi- 
mately 40%)  than  the  electric  compressor  refrigerator. 

The  thermal  compressor  refrigerator  contains  a thermal  compressor,  an 
expander  assembly,  a power  conditioning  equipment  module  for  each  assembly, 
and  Interconnecting  piping.  The  thermal  compressor  Is  as  described  In  the 
preceding  section.  The  expander  assembly  Is  similar  to  the  one  being  built 
In  a current  program  (see  Reference  3).  It  consists  of  a central,  cylindri- 
cal vacuum  shell  which  contains  the  expansion  engines,  heat  exchangers,  cold 
piping,  and  auxiliaries.  A rotary  reciprocating  drive  assembly,  which  con- 
trols the  motions  of  the  pistons.  Is  mounted  on  each  end  of  the  vacuum  shell. 
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with  Thermal 


Refrigerator 


with  Electric 
Compressor 


6 DU  x 49 


14  DIA  x 51 


TABLE  Xll 
SYSTEM  COMPARISONS 


Refrigerator 


Input  Power,  Watts 


Electric  Power 


Compressor 

Expander 


TOTAL  Electric 


Thermal  Power 


Size,  Inches 


Compressor  Assembly 


Expander  Assembly 


Compressor  PCE 


Expander  PCE 


70 

70 

20 

15 

175 


72 

20 

92 

1,064 


12  DIA  x 60 
14  DIA  x 51 
9 x 4 x 39 
9 x 4 x 19 


Weight,  Pounds 

Compressor  Assembly 

155 

Expander  Assembly 

70 

Compressor  PCE 

15 

Expander  PCE 

15 

TOTAL 

255 

^The  sizes  given  are  overall  envelope  dimensions. 


Each  drive  assembly  contains  a rotary  motor,  a gas  spring,  a linear 
actuator,  and  electrical  sensors  for  the  control  system.  From  the  thermo- 
dynamic flowsheet  (Figure  5) , It  can  be  seen  that  the  power  extracted  by  the 
engines  Is  17.5  watts  for  the  first  stage,  4.93  watts  for  the  second  stage, 
and  .790  watts  for  the  third  stage.  It  will  probably  be  desirable  to  split 
the  first  stage  load  between  two,  Identical  8.75  watt  engines  In  order  to 
make  the  engine  capacities  for  the  three  stages  more  nearly  the  same.  This 
will  ease  the  problem  of  selecting  a speed  and  stroke  length  which  will 
result  in  reasonable  dimensions  for  all  three  stages.  For  this  design,  one 
first  stage  engine  and  the  second  stage  engine  are  driven  by  one  drive 
assembly,  and  the  other  first  stage  engine  and  the  third  stage  engine  is 
driven  by  the  other  drive  assembly.  The  two  engines  on  the  same  drive 
shaft  are  mounted  in-line.  The  first  stage  engine  is  driven  by  a shaft  from 
the  drive  assembly,  while  the  colder  engine  is  driven  by  a drive  shaft 
which  penetrates  the  head  of  the  first  stage  engine.  The  heat  exchangers 
and  piping  surround  the  engines  within  the  vacuum  shell.  Cold  gas  from 
the  discharge  of  each  engine  is  piped  in  a vacuum  Insulated  line  from  the 
expander  assembly  to  the  item  to  be  cooled. 

The  power  conditioning  equipment  for  each  assembly  is  Identical  in 
concept  to  the  equipment  presently  being  developed  (see  Reference  3).  It 
consists  of  solid  state  power  and  logic  components  mounted  on  a number  of 
printed  circuit  cards  and  contained  within  an  aluminum  enclosure.  Since 
the  thermal  compressor  contains  two  less  linear  actuators  than  the  electric 
compressor,  there  are  two  less  linear  actuator  Inverters,  with  a corres- 
ponding reduction  in  components.  Additionally,  since  the  electrical  power 
levels  in  the  thermal  compressor  are  less  than  the  electric  compressor,  the 
power  handling  portions  of  the  system  (power  supplies  and  Inverters)  con- 
tain fewer,  more  lightly  stressed  components.  Both  of  these  factors  re- 
sult in  a somewhat  higher  reliability  for  the  thermal  compressor  electronics, 
The  power  conditioning  equipment  for  the  expander  assembly  is  nearly  identi- 
cal to  that  presently  under  development,  except  that  it  handles  less  power 
than  this  equipment. 


The  electric  compressor  refrigerator  also  contains  a compressor 
assembly,  an  expander  assembly,  a power  conditioning  equipment  module  for 
each  assembly,  and  Interconnecting  piping.  The  expander  assembly  and  its 
power  conditioning  equipment  are  Identical  to  their  counterparts  In  the 
thermal  compressor  refrigerator.  The  compressor  assembly  and  Its  power 
conditioning  equipment  are  similar  to,  but  smaller  than,  equipment  presently 
under  development  (see  Reference  3). 


The  system  comparisons  In  Table  XII  did  not  Include  an  allowance  for 
the  weight  or  power  of  the  heat  rejection  system  (liquid  pumps,  radiator, 
and  piping) . Our  experience  with  similar  spaceborne  systems  has  been  that 
the  pumping  power  is  in  the  range  of  5-10%  of  the  heat  rejected  (i.e.,  the 
total  Input  power  to  the  refrigerator),  and  that  the  weight  of  the  heat 
rejection  system  Is  50  lb  per  kilowatt  of  heat  rejected,  if  the  heat  is 
rejected  at  300  K. 
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VII.  CONCLUSIONS 


Based  on  the  analysis  described  in  this  report,  it  is  concluded  that 
the  thermal  compressor  is  a viable  alternative  to  an  electric  compressor 
in  spaceborae,  rotary  reciprocating,  cryogenic  refrigeration  systems.  A 
thermal  compressor  is  a precision  machine,  but  no  extraordinary  technical 
problems  are  foreseen,  either  in  analyzing  machine  performance  or  in  de- 
signing and  fabricating  it.  A thermal  compressor  employs  existing  technology 
throughout,  and  the  machine  evaluated  in  this  study  is  of  a general  size 
which  is  comparable  to  existing  devices. 

The  thermal  compressor  has  a number  of  characteristics  which  are  unique. 
The  most  evident  is  the  fact  that  it  uses  heat  as  the  driving  source.  For 
an  ideal  thermal  compressor,  i.e.,  one  with  no  losses,  the  thermal  input 
power  is  the  isothermal  work  of  compression  multiplied  by  the  efficiency  of 
a reversible  heat  engine  operating  between  the  heat  source  temperature  and 
the  heat  sink  temperature • In  a practical  device , the  thermal  input  power 
is  in  the  vicinity  of  3-4  times  this  power. 

The  volumetric  efficiency  of  a compression  stage  is  a function  of 
the  compression  ratio — with  high  volumetric  efficiency  resulting  from  low 
compression  ratios  and  vise-versa.  For  a single  stage,  the  maximum  volu- 
metric efficiency  is  approximately  60%,  and  the  maximum  pressure  ratio  is 
approximately  2.  In  a practical  compression  stage,  the  volumetric  effi- 
ciency will  be  in  the  range  of  25-35%,  and  the  pressure  ratio  will  lie  in 
the  range  of  1.4  to  1.6.  Thus,  the  thermal  compressor  is  a low  pressure 
ratio,  low  volumetric  efficiency  machine  compared  to  a positive  displacement 
compressor.  Consequently,  several  stages  of  compression  must  be  used  to 
achieve  the  pressure  ratios  required  by  cryogenic  refrigerators  utilizing 
the  reversed  Brayton  cycle. 

Some  electric  power  must  be  supplied  to  the  thermal  compressor  to 
reciprocate  the  displacer,  and  there  is  a trade-off  between  electrical 
power  and  thermal  power.  The  electrical  power  may  be  reduced  to  very  low 
values  if  one  is  willing  to  provide  a relatively  high  amount  of  thermal 
power.  The  electrical  drive  power  is  greater  for  high  speed  machines 
than  for  low  speed  machines.  For  spaceborne  systems,  in  which  one  desires 
to  minimize  electric  drive  power,  reasonable  designs  will  result  from 
reciprocating  speeds  in  the  range  of  300-900  CPM. 

In  comparing  a refrigeration  system  employing  a thermal  compressor 
with  one  employing  an  electric  compressor,  one  finds  that  the  weight  and 
size  of  the  two  systems  are  comparable,  and  that  the  major  differences  are 
in  the  input  power.  Using  the  electrical  power  requirement  of  the  system 
with  an  electric  compressor  as  a baseline,  the  system  with  a thermal  com- 
pressor will  require  15-25%  as  much  electrical  power  and  2-2  1/2  times  as 
much  thermal  power.  The  reliability  of  the  system  with  the  thermal  com- 
pressor appears  to  be  somewhat  higher  than  that  of  a system  with  the 
electrical  compressor,  due  to  a reduction  in  the  component  count  and  elec- 
trical stress  levels  in  the  power  conditioning  equipment. 
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Since  the  thermal  compressor  utilizes  state-of-the-art  components, 
and  since  It  offers  the  opportunity  for  substantial  reductions  In  elec- 
trical Input  power.  It  Is  an  attractive,  viable  alternative  to  an  electric 
compressor  in  spaceborne,  rotary  reciprocating  refrigeration  systems.  It 
Is  recommended  that  thermal  compressor  refrigerators  be  considered  along 
with  electric  compressor  refrigerators  In  spacecraft  system  studies  In  order 
to  determine  the  relative  merits,  from  the  total  spacecraft  system  point 
of  view,  of  the  two  refrigeration  systems.  It  Is  further  recoonended  that 
an  engineering  test  model  of  at  least  one  stage  of  a thermal  compressor  be 
designed,  built,  and  tested  If  system  studies  Indicate  that  the  thermal 
compressor  has  system  advantages. 
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APPENDIX  I 


■PRECXDHO  PA as  BLAUK-NOT  1. 


EQUATIONS  DESCRIBING  THERMODYNAMIC  PROCESSES 
A.  CYLINDER  PROCESSES 

The  mass  flow  throughput  for  a cylinder  per  cycle,  MCY  is  given  by 


MCY  - MC  - MA 

We  define  the  reduced  void  fraction  by 

rvf  = X2V  . TC  V3V  TC  V4V 

V3D  TH  V3D  T4  * V3D 

the  volumetric  efficiency  by 


«v  = m : l • tc 
PS  • V3D 


(Al) 


(A2) 


(A3) 


and  note  that 


V2D 

V3D 


A2 

A3 


(A4) 


Then,  by  expressing  MC  and  MA  In  terms  of  the  pressures,  volumes  and 
temperatures  of  each  element  In  the  system,  (Al)  through  (A4)  can  be 
combined  to  yield 


NV 


(xi  - RVF)  - 1 (I  - 


(A5) 


which  relates  the  mass  flow  throughput  per  cycle  to  the  pressure  ratio* 
MCY  is  related  to  the  total  mass  flow  rate  by 

W 


MCY 


N • NC 


(A6) 


Thus,  if  W,  N,  NC,  PD,  PS,  TC,  TH,  RVF  and  ||  are  specified,  (A3)  and 
(A5)  can  be  used  to  determine  V3D. 

Now 


V3Dn 

* ir  * SDR 


A3 


f“2 


(A7) 

(A8) 
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S - SDR  • D 


(A9) 


-(a) 


(A10) 


Thus,  if  SDR  is  specified,  the  basic  dimensions  of  the  cylinder  can  be 
determined  from  V3D  and  (A7)  through  (A10). 


The  displacer  positions  at  point  b,  where  the  suction  valve  opens, 
and  at  point  d,  where  the  discharge  valve  opens,  are  obtained  from  mass 
balances  between  points  and  b^  and  c^  and  <1,  respectively,  with  the 
result 


[If  - 1)(Ht8VF) 


A3  #A2  TCl 
V3D  \A3  " THJ 


(All) 


XD  - 


i(fl  + RVF)-(l  + RVF) 


A3  /A2  TC\ 
V3D  \A3  " TH/ 


(A12) 


The  relationships  between  displacer  motion,  regenerator  flow  rate, 
inlet  and  outlet  flow  rates  and  pressure  depend  on  the  portion  of  the 
cycle  that  is  occurring,  as  noted  below. 


Re-expansion,  a-b: 


Cylinder  suction  and  discharge  flows  are  zero.  Regenerator  flow 


is  obtained  from 


WR  - M2 


where  the  dot  indicates  a time  derivative.  By  substitution  we  have: 


(P  • V2  + V2  • £) 


(A13) 


P can  be  related  to  displacer  position  by  a mass  balance  on  the  system 
between  point  a,  and  a point  between  a and  b where  the  displacer  is  at 
position  X,  to  yield 


- 


T 


n 


PD  • V3D  /TC 


P - 


TC 


(1 + RVF) 


x (A2  _ A3\  V3D  /TC 
X \TC  TH  J TC  \TH  RVF/ 


(A1A) 


P can  be  obtained  from  (A14),  while  V2  and  V2  can  be  directly  related 
to  the  displacer  position  and  velocity,  so  that  (A13)  can  be  written  as 


PD  • A2  • XDOT 


WR  - 


V2V 

l1  ~ jgr ; U 


R • TC  • (ex  + 1)‘ 


(A15) 


where 


I 1 


(A16) 


XDOT  - ttN  • S • cos  |2ttN  • 6 


(A17) 

(A18) 


6 Is  a time  parameter  chosen  so  that  X ■ 0 at  6*0. 


Suction,  b-c: 


The  pressure  Is  constant  at  PS.  Regenerator  flow  rate  Is  obtained 


from 


WR  - - M3 


which  can  be  written  as 


WR  - 


PS  • A3  • XDOT 


R • TH 

The  suction  flow  rate  Is  obtained  from 

WS  - M2  + M3 

which  can  be  expressed  as 

PS  • A3  • XDOT 


(A19) 


WS  - 


R • TC 


(A2  TC\ 
A3  “ TH/ 


(A20) 


L 


Compression,  c-d: 

Cylinder  suction  and  discharge  flow  rates  are  zero.  P can  be 
related  to  displacer  position  via  a mass  balance  on  the  system  between 
point  c^  and  a point  between  c^  and  <1  where  the  displacer  is  at  position 
X to  yield 


PS  • V3D 
TC 


fef  + RVF) 


(A21) 


Regenerator  flow  rate  is  derived  from  Equation  (A13),  as  in  the  re- 
expansion process,  and  Equation  (A21),  with  the  result 


WR 


PS  • A2  • XDOT  (l  - ^ fe | + RVf| 

(th  + RVF) 


(6X  + 1)‘  r . TC 


(A22) 


Discharge,  d-a; 

The  pressure  is  constant  at  PD.  Expressions  for  the  regenerator 
and  discharge  flow  rates  are  derived  from  the  same  basic  equations  used 
for  the  suction  process  yielding 


WR 


PD  • A3  • XDOT 
R ♦ TH 


(A23) 


and 


WD 


PD  • A3  ♦ XDOT 
R • TC 


(A2  Tc\ 

A3  ~ TH J 


(A24) 


B. REGENERATOR 


from 


LR,  VFR  and  DH  will  be  specified.  Pressure  drop  is  determined 


Thus,  thermal  performance  is  based  on  mean  flow  rate  through  the 
regenerator. 


C.  HEAT  EXCHANGERS 


The  approach  is  to  specify  a configuration,  the  effectiveness,  E, 
allot  to  the  exchanger  a certain  void  volume,  then  determine  its  pres- 
sure drop  as  a function  of  stroke. 


The  approach  AT  in  the  exchanger,  TWGO,  is  given  by 
TWGO 


2 

WRBAR 


and 


NTU  - In 


- QINBAR  (l  - e\ 
IBAR  • CP  \ E / 


For  an  annular  design  (Figure  11),  we  specify  BAN,  then 

NRBAR  - 2 * WRBAR 

NRBAR  BAN  • U 

nst  - hi * 

NRBAR  • NPR 

based  on  fully  developed  laminar  flow. 

Now,  by  definition 

NTU  - NST  • 


(Cl) 

(C2) 

(C3) 

(C4) 


(C5) 


where 


for  an  annulus.  Also 


DH  - 2C 


VHX  - BAN  • C • L 

(C5)  through  (C7)  determine  C,  L and  DH. 

Now 

c - h 

u tniv 


(C6) 

(C7) 


(C8) 
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so  that 


DPHX  - 


96  L • U • V • G 


DH‘  • 2g 

in  which  V and  G are  functions  of  the  displacer  position. 


(c: 


For  a multipassage  design  (Figure  11),  we  specify  LP.  TWGO  and 
NTU  are  determined  by  (Cl)  and  (C2).  Then 


NTU  - NST  • 


4 LP 


DH 


DH 

NST  - 


BP 

2.89 


NRBAR  - 


NRBAR  • NPR 
WRBAR 


(CJ 

(C] 

(C] 


NP  • BP  • U 


(C] 


VHX  - NP  • BP  • LP 
WR 


G - 


NP  • BP 


F 


14.2 

NR 


(Cl 

(Cl 

(Cl 


so  that 

Dm  . 56-8  L • U • V • C (c 

DH2  • 2g 

Again,  V and  G are  functions  of  displacer  position. 

The  viscosity,  U,  In  the  above  equations  Is  based  on  TC  for  the 
ambient  end  heat  exchanger  and  TH  for  the  hot  end  heat  exchanger* 


D.  VOID  VOLUMES 


Void  volume  relationships  include 


VCL2  + VCHX  + VCFP 


WREG  + VP  DC 


LCFP 


AHFP 


LHFP 


ACFP  and  AHFP  are  determined  from  the  same  equation,  i.e 


WRBAR 


ACFP  or  AHFP 


where 


E.  HOT  END  POWER  INPUT 


, is  obtained  from 


which  can  be  shown  to  be 


It  is  interesting  to  determine  the  ratio  which  is 


SI  . 
w 


m*  th  - tc] 


R • TC  • In 


W] 


The  term  in  the  brackets  is  the  isothermal  work  of  compression  between 
A2 

PS  and  PD.  When  - 1,  the  other  term  is  the  reciprocal  of  the  Carnot 
efficiency  of  an  engine  operating  between  TH  and  TC. 

Thus,  the  ideal  heat  input  to  a thermal  compressor  is  the  same  as 
that  to  a reversible  engine  driving  an  Isothermal  compressor.  In 
principle,  then,  the  thermal  compressor  is  a reversible  device.  In 
practice,  a number  of  thermal  loss  mechanisms  prevent  a close  approach 
to  ideal  performance.  These  losses  are  evaluated  by  the  relations  that 
follow. 

Shuttle  Heat  Transfer 


* • KG  • D • S 


QSH  - 5— C-T 


(TH  - TC) 


Displacer  Conduction 


QDIS  - n * KD1LHC  P * ~ " TC) 


Cylinder  Wall  Conduction 


QCW  - n * KCW.iD  * T-~  (TH  - TC) 


Regenerator  Shell  and  Matrix  Conduction 


QRS  - W ^ ~lRDRS  * "R~  (TH  - TC) 


QRM  - KRMl^  — (TH  - TC) 


Insulation  Conduction 


KINS  • AINS 


(TH  - TC) 
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Pumping  Losses 

The  average  mass  flow  rate  into  (and  out  of)  the  radial  clearance 
between  displacer  and  cylinder  each  half  cycle  is  given  by 


MPDOT 


4tt  « N ■ D « C (PD  - PS)  (LHC  + S] 
R(TH  + TC) 


The  heat  transfer  between  this  gas  flow  and  the  walls  is  characterized 
by  the  NTU's,  given  by 


NTU  - 


1.89  n • D • KG  (LHC  + S] 
C • MPDOT  • CP 


(E10) 


in  which  fully  developed  laminar  flow  is  assumed , 
The  pumping  heat  loss  is  then  given  by 


QPL  - 


Total  Hot  End  Power  Input 


MPDOT  ♦ CP  (TH  - TC] 


2 (NTU  + 1) 


(Ell) 


The  total  heat  input  is  the  summation  of  the  ideal  input  and  the 
losses,  viz., 


QINBAR  - Q3  + QREG  + QSH  + QDIS  + QCW 


+ QRS  + QRM  + QINS  + QPL 


(E12) 


F.  DISPLACER  POWER  INPUT 


The  work  done  on  the  displacer  by  the  gas  (and  vice  versa)  each 
cycle  is  comprised  of  two  parts,  that  due  to  the  pressure  difference 
between  hot  and  ambient  end,  WDP,  and  that  due  to  the  area  difference, 
WDA,  given  by 

WDP  - / A2  (P2  - P3)  DX  (FI 


The  associated  powers  are 


- f P3  (A3  - A2) 


PDP  ■ N • WDP 
PDA  - N • WDA 


and  the  total  displacer  power  is  given  by 


PDISP 


In  (FI)  the  pressure  drop  Is  obtained  from 


(P2  - P3)  - - (DPREG  + DPHHX  + DPCHX  + DPHFP  + DPCFP)  vFb) 

In  which  all  the  quantities  are  evaluated  as  functions  of  the  displacer 
position,  X. 

We  also  note  that  (F2)  is  equivalent  to 


Thus,  as  one  would  expect,  the  work  done  on  the  area  (A3 
#A3  - A2\ 

fraction  I 7-5 1 of  the  p-V  work  done  on  A3. 


““ 
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APPENDIX  II 

DESCRIPTION  OF  THERMAL  COMPRESSOR  COMPUTER  PROGRAM 


A.  GENERAL 

4«r  “aly8i8  of  the  behavlor  of  a thermal  compressor  necessitates  the 
specification  of  a large  number  of  parameters  prior  to  making  any  calcu- 
lations. The  calculations  themselves  are  straightforward;  however,  they 
are  tedious  and  would  be  time-consuming  to  perform  by  hand.  We  have  opted 
for  a simple  computer  program  to  do  these  calculations:  optimization  of 

the  machine  Is  then  a matter  of  judgment,  requiring  the  exploration  of 
several  of  the  more  critical  Input  parameters  over  their  appropriate  ranges. 

The  program  is  a simple,  single-stream  computational  sequence  as 
shown  in  Figure  21.  It  consists  of  a main  program  and  two  subroutines: 
one  for  the  regenerator  and  one  for  the  heat  exchanger  designs.  Execution 
proceeds  as  follows:  the  forty-one  input  parameters  are  read  in  from  card 

input  and  printed  with  their  respective  descriptions  on  a title  page.  This 
page  serves  to  completely  define  the  problem,  and  it  eliminates  the  need 
or  a title  card.  Next,  the  program  calculates  the  geometric  aspects  of 
the  thermal  compressor  (cylinder  diameter,  void  volumes,  port  locations) 
and  also  calculates  the  displacer  position  and  velocity  at  one  hundred 
increments  of  the  cycle.  At  the  same  time,  the  ambient  end  pressure  and  the 
regenerator  flow  rate  are  evaluated  at  each  of  these  100  displacer  positions. 

Once  t le  regenerator  flow  rates  are  defined,  the  regenerator  design 
subroutine  Is  called  and  a thermal  design  Is  obtained  using  mean  flow 
properties.  After  this  step,  the  pressure  drop  (based  on  instantaneous 
flow)  across  the  regenerator  is  computed  at  each  of  the  100  cycle  increments. 
In  a similar  fashion,  the  heat  exchanger  subroutine  is  called  to  design  the 
hot  and  ambient  end  heat  exchangers. 

Returning  to  the  main  program,  we  proceed  to  evaluate  thermal  losses 
for  the  machine  using  the  calculated  geometry  and  the  input  values  of  thermal 
conductivities,  wall  thicknesses,  etc.  The  last  event  in  the  computational 
sequence  is  a stepwise  trapezoidal  integration  of  the  pressure  drops  in  the 
regenerator  and  the  two  heat  exchangers.  The  total  loss  through  these  three 
elements  must  be  offset  by  a shaft  power  input  to  the  machine.  This  power 

input  is  simply  a product  of  the  speed  and  the  integral  of  the  total  pressure 
drop  over  a single  cycle. 

The  nomenclature  used  in  the  program  is  defined  in  the  list  of  abbre- 
viations and  symbols. 


B.  INPUT  DATA 

Input  data  for  the  program  consist  of  41  parameters  on  6 cards.  Defi- 
nition of  these  variables  and  the  required  format  on  the  input  computer  cards 


..  . 


are  shown  In  Table  XIII  In  the  order  In  which  the  cards  are  read.  Units 
for  all  variables  are  listed,  along  with  the  proper  location  of  each 


parameter  on  the  data  cards.  There  Is  no  title  card  Included  In  the  In- 
put since  virtually  all  of  the  Input  parameters  are  listed  In  the  output, 
thereby  completely  defining  the  problem.  Some  of  the  Input  parameters  are 
discussed  in  the  following  paragraphs. 


1.  Void  Fractions 

All  void  fractions  are  defined  as  the  respective  void  volume  relative 
to  the  hot  end  swept  volume. 

The  hot  and  ambient  end  void  fractions  account  for  the  clearance  volume 
at  cylinder's  end,  for  the  heat  exchangers,  and  for  the  gas  flow  passages. 
Porting  void  fraction  Is  the  dead  volume  of  the  ambient  end  ports.  It  Is 
also  Included  In  the  ambient  end  void  fraction,  but  It  must  also  be  listed 
as  a separate  entry. 

Intermediate  void  fraction  Includes  the  regenerator  void  volume  and 
the  annular  clearance  between  the  displacer  and  the  cylinder. 


2.  Displacer  Area  Ratio 

On  Card  2 the  ratio  of  ambient  displacer  area  to  hot  displacer  area 
must  be  specified.  The  deviation  of  this  value  from  unity  Is  In  recogni- 
tion of  the  presence  of  a shaft  on  the  ambient  end  of  the  displacer.  Hence, 
this  area  ratio  Is  always  less  than  one. 


3 . Cylinder  Length 

The  cylinder  length  Is  the  length  of  the  cylinder  over  which  the  hot 
end-to-amblent  temperature  difference  occurs. 

4.  Regenerator  Configuration 

Column  4 of  the  third  Input  card  sets  a flag  which  determines  the 
regenerator  configuration.  The  two  choices  available  are  (1)  an  annular 
regenerator  which  Is  Integral  with  the  cylinder,  and  (2)  a cylindrical 
regenerator  which  Is  separate  from  the  cylinder.  In  the  latter  case, 
the  volume  of  flow  passages  to  the  regenerator  must  be  Included  as  part 
of  the  hot  end  and  ambient  void  volumes.  Neither  of  the  regenerator 
configurations  has  a unique  matrix  material;  rather,  one  of  three  matrices 
(spheres,  screens,  and  round  tubes)  may  be  selected  by  the  appropriate 
flag  in  Column  8 of  the  same  card. 
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TABLE  XIII 


INPUT  DATA  TO  PROCRAM 


Representative 

Values 


CARD  1:  FORMAT  (10F8.0) 


TH 

TC 

PD 

PS 

N 

NC 

W 

VFH 

VFC 

VFI 


High  temperature  source  (K) 

Low  temperature  sink  (K) 

Absolute  discharge  pressure  (atm)  j 

Absolute  suction  pressure  (atm)  1 

Speed  (Hz) 

Number  of  cylinders 

Total  mass  flow  rate,  all  cylinders  (g/sec) 
Hot  end  void  fraction 
Ambient  end  void  fraction 
Intermediate  void  fraction 


800  - 
300  - 


1,000 

350 


1.25 


PD 

" ps  < 2.°; 


.5  < PS  < 4 


5 

1 

0.2 

.05 

.05 

.3 


20 

2 

1.0 

.25 

.25 

1.0 


CARD  2:  FORMAT  (10F8.0) 


TCW 

TD 

TINS 

XCL3 

XCL2 

SDR 

A2A3 


LHC 

VFPORT 


Cylinder  wall  thickness  (cm) 

Displacer  wall  thickness  (cm) 

Insulation  thickness,  hot  end  (cm) 

End  clearance,  hot  end  (cm) 

End  clearance,  ambient  end  (cm)  — • 
Stroke-to-dlameter  ratio 

Ratio  of  ambient  displacer  area  to  hot  displacer 
area;  this  ratio  is  less  them  one 

Radial  clearance  between  displacer  and  cylinder 
wall  (cm) 

Cylinder  length  between  hot  and  ambient  ends  (cm) 
Porting  void  fraction 


.025 

.05 

3 

.05 

.05 

.1 


.050 

.10 

10 

.10 

.10 

.5 


.9  - .95 


.03 

5 

.01 


.10 

20 

.05 


CARD  3:  FORMAT  (214,  9F8.0) 


IREG 

IMTX 


In  Column  4: 


In  Column  8: 


DH 

LR 

VFR 

TRS 


1 - regenerator  is  Integral  with  cylinder 

2 - regenerator  is  separate  from  cylinder 

1 ■ spheres 

2 ■ screens 

3 ■ tubes 

Hydraulic  diameter  of  regenerator  matrix  (cm) 
Regenerator  length  (cm) 

Regenerator  void  fraction 
Regenerator  shell  thickness  (cm) 


.01 

1 

.02 

.025 


.10 

10 

.10 

.050 
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TABLE  XIII 


CARD  4: 

IHHX 

LPH 

EHHX 

LHFP 

DYNNH 

CARD  5: 

ICHX 

LPC 

ECliX 

LCFP 

DYNHC 

CARD  6: 

KCW 

KDIS 

KRS 

KINS 

KRMM 


(Concluded) 


FORMAT  (14,  9F8.0)  HOT  END  HEAT  EXCHANGER 

In  Column  4:  1 ■ annular  heat  exchanger 

2 “ multipassage  heat  exchanger 
Passage  length  In  multipassage  design  (cm) ; leave 
blank  In  annular  design 
Effectiveness 

Length  of  fluid  passage  between  hot  end  and 
heat  exchanger  (cm) 

Dynamic  head  In  fluid  passage  (atm) 

FORMAT  (14,  9F8.0)  AMBIENT  END  HEAT  EXCHANGER 

In  Column  4:  1 » annular  heat  exchanger 

2 - multipassage  heat  exchanger 
Passage  length  In  multipassage  design  (cm);  leave 
blank  In  annular  design 
Effectiveness 

Length  of  fluid  passage  between  ambient  end  and 
heat  exchanger  (cm) 

Dynamic  head  In  fluid  passage  (atm) 

FORMAT  (10F8.0) 

Thermal  conductivity  of  cylinder  wall 
(watt/cm-K) 

Thermal  conductivity  of  displacer  wall 
(watt/cm-K) 

Thermal  conductivity  of  regenerator  wall 
(watt/cm-K) 

Thermal  conductivity  of  insulation  on  hot  end 
(watt  cm-K) 

Thermal  conductivity  of  regenerator  matrix 
material  (watt/cm-K) 


Representative 

Values 


1-10 
.8  - .9 

0-15 
.01  - .1 


1-5 
.8  - .9 

0-5 
.01  - .1 

.1  - .2 
.1  - .2 
.1  - .2 
0001  - .001 
.1  - .2 


5.  Regenerator  Void  Fraction 


Porosity  of  the  regenerator  matrix  material  must  be  taken  Into  account 
by  specifying  the  matrix  void  fraction  in  Columns  33-40  of  Card  4. 


6.  Choice  of  Heat  Exchanger 

If  a multipassage  heat  exchanger  is  selected,  the  length  of  a single 
passage  must  be  specified  in  Columns  5-12  of  Card  4 or  5.  If  an  annular 
configuration  is  selected,  there  need  not  be  an  entry  for  length.  The 
program  will  set  the  perimeter  of  the  annulus  equal  to  that  of  the  cylinder 
for  the  annular  heat  exchanger. 


7 . Fluid  Passages 

The  fourth  and  fifth  data  entries  on  Cards  4 and  5 specify  the  length 
of  the  fluid  passages  to  the  heat  exchangers  and  the  dynamic  head  in  those 
passages.  Together  with  the  calculated  regenerator  flow,  these  parameters 
will  determine  the  cross-soctlonal  area  of  the  fluid  passages.  The  volumes 
of  the  passages  are  then  subtracted  from  the  total  void  volume  to  yield 
the  volume  available  for  the  heat  exchangers. 


C.  HELIUM  PROPERTIES 

We  have  modeled  the  working  fluid  as  an  ideal  gas  since,  in  the 
range  considered  (2-4  atmospheres,  300-900  K) , the  gas  is  indeed  ideal, 
and  furthermore,  has  a constant  specific  heat  and  Prandl  number.  Fitting 
the  thermal  conductivity  exponentially  in  temperature  then,  we  obtained 
dose  agreement  with  published  NBS  data  (Reference  10) . The  viscosity  also 
has  a functional  dependence  on  temperature.  The  constant  of  proportionality 
was  derived  from  the  specific  heat  and  Prandl  number. 


D.  REGENERATOR  DESIGN 


1.  Heat  Transfer  and  Friction  Data 

Three  matrix  configurations  are  incorporated  in  regenerator  sub- 
routine REGDES:  spheres,  screens,  and  round  tubes.  Wherever  possible, 

the  analytical  correlations  provided  by  Kays  and  London  (Reference  11) 
are  used.  In  cases  where  correlations  were  not  available,  the  graphical 
or  tabular  data  in  Kays  and  London  were  curve-fitted  in  Chebyshev  polynomials. 

In  the  case  of  tubes,  virtually  any  porosity  may  be  specified,  making 
the  algorithm  quite  general  in  that  respect.  On  the  other  hand,  the  corre- 
lation for  spheres  is  valid  only  for  porosities  in  the  0.37-0.39  range 
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(the  maximum  theoretical  packing  density  yields  a porosity  of  0.37);  and 
the  correlation  for  screens  is  valid  only  for  the  minimum  porosity  (0.65) 
for  which  Kays  and  London  provide  data. 


Shell  Conduction 


Conduction  through  the  regenerator  shell  is  straightforward.  Two 
regenerator  configurations  are  allowed:  (1)  an  annular  regenerator  which 

is  integral  with  the  cylinder  wall,  and  (2)  a separate  external  regenerator 
connected  to  the  heat  exchangers  via  piping.  For  the  annular  regenerator, 
two  shells  are  required  since  we  allow  for  insulation  between  the  inside 
of  the  regenerator  and  the  cylinder  wall.  This  arrangement  permits  us 
to  choose  a regenerator  length  which  is  shorter  than  the  cylinder  length 
and  still  maintain  separate  temperature  profiles  in  each.  The  diameter 
of  this  inside  shell  is  set  to  be  1.0  inch  larger  than  the  cylinder  diameter. 
Thermal  paths  through  both  shells  are  accounted  for  in  the  conduction  loss. 
Should  the  regenerator  and  cylinder  lengths  be  chosen  equal,  the  insulation 
is  not  required  and  a single  regenerator  shell  would  suffice. 


. Regenerator  Matrix  Conduction 


We  have  invoked  a rigorous  analysis  of  matrix  conduction  effects  only 
in  the  case  of  round  tubes.  Given  a hydraulic  diameter  and  regenerator 
void  fraction,  the  tube  diameter  and  wall  thickness  may  be  calculated  if 
the  tubes  are  arranged  in  a hexagonal  packing.  The  actual  conduction  loss 
may  be  calculated  from  this  result  given  the  frontal  area  and  length  of  the 
regenerator. 


In  the  case  of  spheres,  we  based  the  loss  formulation  on  n cubic 
packing,  i.e.,  a packing  with  the  spheres  centered  at  the  corners  of  a 
cubic  lattice.  Inspection  of  a single  stack  of  spheres  reveals  that  the 
thermal  path  is  simply  a series  of  contact  resistances,  assuming  that 
the  thermal  resistance  through  the  sphere  itself  is  negligible.  Ws  have 
assumed  a mechanical  loading  of  0.068  atm  (1.0  psi)  on  the  spheres  as  a 
deterrent  to  matrix  movement.  For  the  total  matrix,  then,  multiple  con- 
ductive paths  are  present,  the  number  being  a function  of  sphere  diameter 
and  regenerator  frontal  area.  There  is  an  inconsistency  in  this  conduction 
model  since  our  heat  transfer  and  friction  data  is  based  on  a matrix 
porosity  of  0.37-0.39,  which  would  require  a "hexagonal  close  packed" 
lattice.  The  cubic  packing  yields  a porosity  of  0.48.  Our  principal  justi- 
fication for  this  approach  to  calculating  matrix  conduction  is  its  simplicity. 
It  spite  of  the  inconsistency,  it  is  at  least  a tacit  recognition  of  the 
existence  of  matrix  conduction. 


Screens  are  no  easier  to  analyze  than  spheres.  Hence,  our  rationale 
in  evaluating  solid  conduction  through  a matrix  of  screens  follows  very 
closely  the  rationale  applied  to  spheres,  viz.,  in  the  interests  of  simpli- 
city, try  a formulation  which  will  be  cognizant  of  the  phenomenon  and 
which  will  yield  some  useful  results.  Accordingly,  we  have  modeled  the 
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screen  as  a "crossed-rod"  lattice  with  the  rods  In  each  layer  oriented  at 
90°  to  those  In  the  adjacent  layer.  The  mesh  size  and  wire  diameter  of 
this  configuration  may  be  derived  solely  from  a matrix  hydraulic  diameter 
If  one  assumes  that:  (1)  each  layer  Is  two  wire  diameters  thick,  and 

(2)  the  wire  spacing  Is  the  same  for  all  layers  (Reference  12) • Correla- 
tions of  wire  diameter  and  wire  density  as  functions  of  hydraulic  diameter 
are  used  for  this  calculation.  A series  of  contact  resistances  at  the 
Intersection  of  two  wires  then  form  the  thermal  path  through  the  matrix. 
Given  the  regenerator  length,  we  can  determine  the  number  of  resistances 
i-i.  series;  the  frontal  area  will  then  yield  the  number  In  parallel. 


E.  HEAT  EXCHANGER  DESIGN 

In  subroutine  HXDES,  a heat  exchanger  for  either  the  hot  or  aublent 
end  Is  designed.  Vo  configurations  are  permitted:  (1)  annular,  in  which 

the  heat  exchanger  consists  of  an  uninterrupted  annular  space  In  the  cylinder 
wall;  and  (2)  multipassage,  In  which  the  annular  space  Is  broken  up  Into 
a number  of  square  passages  around  the  perimeter  of  the  cylinder.  In  each 
of  these  cases,  there  Is  a provision  for  a "fluid  passage"  to  and  from  the 
heat  exchanger.  The  actual  volume  available  for  the  heat  exchanger  Is  the 
void  volume  for  the  particular  end  (hot  or  ambient)  being  examined,  minus 
the  volume  of  the  end  clearance  In  the  cylinder  and  the  volume  of  fluid 
passage (s) . 

For  an  annular  design,  the  program  sets  the  annulus  diameter  equal  to 

Ithe  cylinder  diameter  plus  1.0  Inch,  and  calculates  a perimeter  from  that. 

As  In  the  regenerator,  thermal  design  Is  based  on  mean  flow  properties, 
and  the  pressure  drop  Is  based  on  Instantaneous  flow.  In  the  printout, 
the  perimeter,  length,  and  width  of  the  annulus  are  given  along  with  the 
thermal  design  parameters. 

The  procedure  for  a multipassage  design  Is  similar  to  that  for  the 
annular  design,  viz.,  thermal  design  Is  based  on  mean  flow  properties, 
pressure  drop  on  Instantaneous.  The  printout  lists  the  passage  length 
(Input  parameter),  passage  width,  number  of  passages,  and  the  thermal 
design  parameters. 


F.  DISPLACER  POWER 

At  the  end  of  the  main  program,  a trapezoidal  Integration  technique 
Is  used  to  evaluate  the  shaft  power  required  to  overcome  the  pressure 
drop  losses  In  the  regenerator,  heat  exchangers,  and  fluid  passages.  It 
should  be  remembered  that  we  have  evaluated  the  Instantaneous  flew  through 
the  regenerator  at  100  points  In  the  cycle.  These  100  points  are  distri- 
buted In  four  regimes  (25  points  to  each) : (1)  reexpansion,  (2)  suction, 

(3)  compression,  and  (A)  discharge.  The  \ 'es  .re  drops  in  the  regenerator 
and  In  each  heat  exchanger  and  fluid  passage  a'.e  calculated  at  each  instan- 
taneous flow  rate  and  then  stored  In  arrays.  From  these  values,  the  pressure 
drop — displacer  area — stroke  increment  product  Is  summed  around  one  cycle 
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to  obtain  the  displacer  shaft  work  required  for  a single  cycle.  Displacer 
power  is  simply  the  product  of  this  work  per  cycle  and  the  speed  of  the 
machine. 

Work  output  may  be  derived  from  the  thermal  compressor  as  a result 
of  the  difference  in  hot  and  ambient  displacer  areas  due  to  the  presence 
of  a shaft  on  the  cold  end.  This  work  output  may  be  calculated  in  a 
straightforward  manner  from  the  ratio  of  displacer  areas  and  the  ideal  hot 
end  power  input.  Subtracting  this  value  from  the  pressure  drop  losses 
yields  the  net  displacer  shaft  power  which  will  have  to  be  furnished  to 
the  machine  by  the  linear  actuators.  The  printout  lists  all  of  the  contribu- 
tions to  the  net  displacer  shaft  power  requirement,  as  well  as  the  flow 
rate  and  pressure  drops  through  the  regenerator,  heat  exchangers,  and 
fluid  passages  at  each  of  the  100  stroke  positions  in  a single  cycle. 


G.  LISTING  AND  SAMPLE  PROBLEMS 

A listing  of  the  program  is  shown  in  Figure  22,  and  two  sample  problems 
are  shown  in  Figures  23  and  24.  The  sample  problems  correspond  to  the  final 
design  point  which  was  described  in  Section  V. 
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ThtS  PROGRAM  OFSIGNS  a ROTARY*RbCiPROCATING  THERMAL  COMPRESSOR , 
THE  PROGRAM  OUTPUT  WILL  LIST  ALL  OF  THE  INPUT  OATA  AS  WELL  AS  THF 
GEOMETRIC  ASPECTS  OF  THE  MACHINE.  ITS  THERMAL. AND  SHAFT  POWER 
REOUIREMFNTS,  AND  THF  LOSSES  ASSOCIATED  WITH  THE  DESIGN, 

**♦  CASE  NO,  77056  ( WP AF B»AF POL ) »** 

***  P,M,  O'FARRELL.  18  JULY  74  *** 


REAL  kcw,kdis,kg,kins,krs,lhc,lcfp,lhfp,lr,lph,lpc»krmm 

REAL  N.NC.NPR 

REAL  NV,MCY,MPD0T,NTU 

DIMENSION  OPHHX(100),OPCHX(100),X1  (10O),,TYP(16),WCHX(12) 

DIMENSION  OPHFP(100),OPCFP(100) 

COMMON/ INPUT/TH, TC.TCR# PS, PO.N 

COmmon/RFG/WREG,  VFR.  V3D.LR,0H,DPREG(  100)  |QREG.QRS.QRM,KRS,TRS. 
IIREG.IMTX.D.LHC.wrmm 
COMMON/SAVE/WR(100) ,P(100) 

DATA  R, CP, NPR. PI/2,079. 5,193. ,667,3,1415926536/ 

DATA  TYP/4HANNU.4HLAR  , 4HSEP A , 4HR ATE , 4HSPHE , 4HRE S , 4H3CRE , 4HENS  , 
14HTUBE.4HS  , 4HANNU, 4HLAR  ,4H  , 4HMUL T , 4HI P AS , 4HS AGF/ 

DATA  WCHX/4HPERI,4HMFTF,4HR(CM,4H)  ,4H  ,4H  , 4HLENG, 4HTH  0, 

1 4HF  PA,4HS3aG,4HE(CM,4H)  / 

READ  THE  INPUT  OATA 

2 RE  AD (2, 500 ) Th, TC ,PD , PS ,N ,NC , w , VFH, VFC . VF I 

RE  AD (2, 500)  TCW,T0,TINS,XCL3,XCL2,S0R,A2A3,C,LHC, VFPORT 
RE AD (2. 501)  1RFG, IMTX.DH.LR.VFR.TRS 
RE  AD (2, 5P2)  !HHX,LPH,£HHX,LHFP,0YNHH 
READ (2, 502)  ICHX.LPC.ECHX.LCFP.DYNHC 
RE AD (2, 500)  KCW,K0!S,KRS,KINS,KRMM 


PRINT  THF  INPUT  DATA 


WRITE (5,700) 
WRITE(5,701 ) 
WRITE(5, 702) 
WRI TE (5, 703) 
WRI TF (5,704) 
WRITE (5, 705) 
WRI TE (5, 706 ) 
wrITE(5,707) 
WRITE (5, 708) 
WR I TE (5, 709 ) 
WRITEf5,7l0) 
WRITE ( 5, 7 1 1 ) 
WRITE (5,712) 
WRITE (5, 713) 
IHXbIHHX 


TM,TYP(2*IREG*1).TYP<2*IREG) 
TC,TCW,TYP(2*IMTX*3) ,TYP(2*IMTX*4) 
PD,TD,DH 
PS, TINS, LR 
n,  VFR 

NC, XCL3.TRS 

W.XCL2 

SDR 

VFH. A2A3 
VFC.C 

VFI,  LHC, VFPORT 


WRITE (5, 714)  KLW,TYP(3*IHX*8),TYP(3*IHX*9) ,TYP(3*IHX410), 
1TYP(3*ICHX*8),TYP(3*ICHX*9),TYP(3*ICHX*10) 

GO  TO  (10,U),THHX 
10  WRITE(5,715)  KDIS 


Figure  22.  LISTING  OF  COMPUTER  PROGRAM 
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GO  TO  12 

It  WRITF(5(716)  KOlS,LPH,(NCHX<6*lCHX.6*IW),IN»l,0>(LPC 
12  WRITE (5,717)  KRS,  EHHX , F.GHX 
WRITE (5, 7 10)  FINS 
WRITE(5#719)  *RMM,LHFP,LCFP 
WRITE(5,723)  DYNHH,DVNHC 
C 

C SIZE  THE  CYLINDER 
C 

RVFpVFC*VFH*TC/TH*VFI*TC/TH 

TCR«TC/TH*RVF 

AT»A2A3*TC/TH 

AR«A2A3*RVF 

NV«A2A3*RVF«TCR*PD/PS 

FG«2,861F»i5*(t5*(TH*TC))**,700 

MCYpW/N/NC 

V30»MCY*R*TC/(PS*NV)/t 101325 

V2V*VFC*V30 

V3V«VFH*V30 

V4V*VFI*V3D 

Da  (4,*V3Q/P1 /SDR) ** ,3333 

SaO*SOR 

A3aPI*D*D/4, 

A2*A3*A?A3 

VCL2aA2*XCL2 

VCL3»A3*XCL3 

VP0CaP!*0*C*(LHC*8) 

VVREGbV4V*VPDC 

C DISPLACER  POSITIONS  AT  OPENING  OF  SUCTION  ANO  OISCHARGE  PORTS 

C 

XBi(PD/PS"l,)*TGR/!A3aAT/V3D) 

X0a((PS/P0)*AR.TCR)/(A3*AT/V30) 

C 

c EVALUATE  FLOW  rate  and  pressure  versus  POSITION 

c 

c REEXPANSION  a*b 

c 

DO  50  1*1*25 
Xal*XB/25( 

XKlIaX 

XD0TaS*PI*N*SQRT(l,«(2(*X/S«ll)**2) 

ZaV3D*TCR/TC 

P(I)aPD*Z/(X*(A2/TC»A3/TH)*Z) 

Z*(PD/R/TCT *TCR*(T2R«A2A3*VFC*AT)*A2*X00T 
Zl«(A3*X*AT/V3D*TCR)**2 
wR(I)a,10132S*Z/Zl 
50  CONTINUE 
C 

C SUCTION,  B-C 
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0090 

0099 

0100 


0101 
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0103 


0105 


0106 

0107 


NR(naP3»A3*XDOT*,10132b/(R*TH) 
P(I)iPS 
60  CONTINUE 


COMPRESSION,  C*D 


00  70  H51,75 
XaS-(I»50)*(S«XD)/25, 

XI  max 

IF(XtCT(S)  Xa0 

X0nTaS*Pl*N*S0RT( l ,* ( 2 , *X/S- 1 , ) * *2) * (» 1 , ) 

P(I)aPS*AR/(X*A3*AT/V3D*TCR) 

Za (PS/R/TC) *AR* (TCR-VFC*A2A3*AT) *A2*XD0T 
ZlaX*A3*AT/V3D*TCR 
mR ( I ) ■ , 101 325*Z/Z 1 
70  CONTjNUt 


DISCHARGE t D*A 


DO  00  Ia76»100 
XaXD*  f 1 , •(1*75)725,) 

X 1 ( I ) a X 

X0OTaS*PI*N*S0RT ( 1 ,*C2,*X/S-1 ,)**?)♦ (»1 , ) 
WHmaPD*A3*X00T*,  101325/(R*TH) 

P(I)aPD 
00  CONTINUE 


DESIGN  THE  REGENERATOR 


WRI TE (5, 724 ) 

HHITE (5, 720)  TYP(2*IREG*l),TVP(2aIREC) 
CALL  REGDES 


CALCULATE  THE  IDEAL  HEAT  INPUT  PER  CYCLE 
Q3a, 101 325*NV*PS*V3D*AL0G (PD/PS) /AT 
03*N*Q3 


SHUTTLE  LOSS 


QSHaKG*PI*D*S*S* (TM..TC  )/ (5*4*C*LMC) 


CONDUCTION  THROUGH  DISPLACER 


ODISaKDIS*PI*D*TD*(TH»TC)/LHC 


CONDUCTION  THROUGH  CYLINDER  HALLS 


QCWaKCMaPl*0*TCH*(TH*TC)/LHC 


conduction  THROUGH  insulation 


AINSaPI*fD*TINS)*( (D*TINS)/4,«S*LR*TIN8) 
QIN3«KINS*AIN0*(TH. TCI /TINS 


Figure  22  (Continued) 
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0108 

0109 

0110 
0111 


0112 


PUMPING  LOSSES 


MPDOTM,»N*PI»0*C*(PD*PS)*(LHC*S)/R/(TH*TC) 
MPOOT  * MPDOT* ( 101 32S 
NTU>1(89»PI»D»(IHC«S)»KG/(C*MP00T«CP) 
OPL«MPDOT*CP*(TH»TC)/C2i«(NTU*1,)) 

total  hot  end  POWER  input 

0INBAR>G3*QRFG*0SH«QDI3«QCW»QRS«QRM4QINS4GPL 


DESIGN  THE  HEAT  EXCHANGERS 
HOT  END  HX 


9113 

0114 


011b 

0118 

0117 

0118 


0119 

0120 
0121 
0122 

0123 

0124 

0125 

0128 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 
0136 

0139 

0140 

0141 


HOT  END  HX 

WRITF (5# 721)  TYP(3*IHX*8),TYP(34IHX*9),TYP(3«IHX410) 
CALL  HXDES(TH,QINBAR,EHHX#V3VtVCL3#IHHX|DPHHX,LPHtLHPP* 
1FP) 


OYNHH,0«OPH 


ambient  end  hx 


IHXxICHX 

WRITF (5, 722)  TYP( 3* IHX*8) , TYP(3* IHX*9) , TYP (3*IHX«1«) 

SUBTRACT  THE  PORTING  VOLU/r  FROM  THE  TOTAL  COLO  VOID  VOLUME 
V2VPbV2V*VFP0RT«V3D 

CALL  HXDFS(TC,0INBAR,ECHX,V2VPlVCL2.ICHX,DPCHX#LPC»LePP«0YNHC,0,DP 
1CFP) 

EVALUATF  D13PLACFR  POwER 
ADIFbA3*A2 

P2P3»»(DPREG( l )*DPHHX ( l )*DPCHX ( l ) ) 

WDPa(6*A2*(P2P3*0l)*(X6/25()*l 101325 
ZiP2P3 

P3»P( 1 )»P2P3 

NOA««I5*ADIF*(PD*P3)*DX*(101325 
Zl  »P3 

DO  150  I *2* 100 

P2P3»»  (OPREG  (I)*DPHHX(1)  *DPCHX  ( I ) *DPHF  P(l)*OPCFP(D) 

OXbXB/25, 

IF ( I ,GT ,25)  DXa(S*XB)/25( 

IF ( I ,GT,50)  DX»»tS»X0)/25( 

IF (l,GT,7b)  0X«*XD/25, 

HDP«W0P*,5*A2* (P2P3*Z) *DX* , 101 325 
P3»P ( I )*P2P3 

WDA«W0A»,5*ADIF*(P3*Zl)*0XB, 101325 

Z>P2P3 

Zl«P3 

150  CONTINUE 
PDPb*N*nDP 
PDAbNbWQA 
PDAb ( 1 ,*A2A3) *Q3 
POISPtPOP-POA 

Figure  22  (Continued) 
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PRINT  THF  COMPRESSOR  GEOMETRY  AND  POKER  REQUIREMENTS 

MR  I TF (5,730)  RVF,Q3,NV,QSM,mCV,Q0IS#V30#QCK 
WRITE (5,731)  V2V,QInS,V3V,QPL,V4V,D,QREG 
WH I TF ( 5 , 73? ) S,QRS, A3,QRM, A2,QINBAR, VCL2 
WRITE (5,733)  VCL3 , VPDC , VVREG, POP, XU, PDA 
WRITE (5,734)  XD,PDISP 

# 

PRINT  THE  PRESSURE  DROPS  AT  FACH  CYCLE  INCREMENT 

WRITE (5,510) 

DO  200  Iwal ,100 

200  WRITE (5, 601)  IW,X1(IW),P(IM), wR( In) ,DPREG( Iw) , DPHHX ( Iw) , OPCHX ( Iw) , 
10PHFP(Ih),DPCFP(IW) 

GO  TD  2 

500  FORMAT! 10F8.0) 

501  FORMAT(2I4,9FB,0) 

502  FORMAT (14, 9F8.0) 

510  FORMAT (/T 10, ' X' ,T22, »P2' ,T34, 'WR1 ,T45, ' DPREG' , T57, 'DPHHX ' ,T69, 'DPC 
1HX«,T82, 'DPHFP»,T94,  'DPCFP'/) 

550  F0RMAT(/2X, 'INPUT  DATA'/) 

600  F0RMAT(/8(1PE12,4)) 

601  FORMAT (2X, I3,8(1P£12,4) ) 

700  F0RMAT(1H1,T45, 'THERMAL  COMPRESSOR!  INPUT  DATA'//) 

701  FORMAT ( T 1 6, ' ••  CYCLE  **',T50,'*«  GEOMETRY  **tfT92,!**  REGENERATOR 
l ** ' /) 

702  FORMAT (T2, 'HOT  END  TEMP ( K ) • , T25, l PE  1 2 1 4 , T39, ' TH ICKNESS (CM)  • , T|4 , A4 
1,A4,T93, 'CONFIGURATION') 

703  F0RMAT(T2, ' AMBIENT  END  TEMP (K ) t , 

1 T25,1PE12,4,T41, 'CYLINOER  WALL»,T66,1 

1PE 12,4, T84,A4,A4,T93, 'MATRIX') 

704  FORMAKT?, 'DISCHARGE  PRESSURE  ( ATM)  ',  T25,  IRE  1 2,4  , T4 1 ,«  0 1 (PLACER  HAL 
1L ' 

1,T65,1PE12.4,T84, 'HYDRAULIC  DI A (CM) ' , T 1 05, 1 PE  1 2 , 4 ) 

705  F0RMAT(T2, 'SUCTION  PRESSURE ( ATM) • , T25, 1 PE  1 2, 4 , T4 1 ,' INSULATION  !, 
1T65,1PE1?,4,TB4, 'LENGTH(CM) ',T105,1PE12,4) 

706  FORMAT (T2, 'SPEFD(HZ) • , T25, 1 PF 1 2, 4 , T39, « AXI AL  CLEARANCE (CM) • ,184, » V 
1010  FRACTION', T105.1PE12, 4) 

707  F0RMAT(T2, 'NUMBER  OF  C YL I NDERS ' , T25, 1 PE  1 2, 4 , T4 1 , « MOT  END ' , T65, 1 PE i 
12,4,TB4, 'SHELL  THICKNESS ( CM ) « , T 1 05, 1 PE  1 2, 4) 

708  FORMAT!!?, • TOTAL  FLOW  RA TE ( G/3EC ) » , T25, 1 PEI 2, 4 , T4 1 , • AMBIENT  END', 
1T65,1PE12.4) 

709  FORMAT ( T2, 'VOID  FRACTIONS ', T39, t STROKE/DI AMETER  RATIO ', T66, IP£1 2,4 
1) 

710  FORMAT ( T6, 'H0T',T25,1PE12,4,T39, 'AREA  RATIO|  AmBIENT/HOT «, T65, 

1 1 PF 1 2,4 ) 

711  F0RMAT(T6, ' AMBIENT (INCL,  PORT)'  , T25, 1 PE  12, 4 , T39 , 1 RA0I AL  CLEARAN 
1CE(CM) i,T65,1PE12,4) 

712  FORMAT (T6, 'INTERMEDIATE • , T25 , 1PE 1 2, 4 , T 39, 'CYLINDFR  LENGTH (CM) « ,T65 
1,1PE12,4/T6, 'PORTING',  T25,1PE12,4) 

713  FORMAT (//T8, •**  THERMAL  CONDUCTIVITY  •*',T46,'*e  HOT  HEAT  EXCHANGE 
1 R *»',T8B,'*»  AMBIENT  HEAT  EXCHANGER  ••'/) 

714  FORMAT (/T2, 'CYLINDFR  WALL(W/CM«K) I , T25, 1 PE  1 2 ,4 , T 39 , 3 A4 , T52, • T YPE ' , 
1T84,3A4,T98, 'TYPE') 

715  F0RMAKT2, 'DISPLACER  WALL  ' , T25, 1 PE  1 2,4 ) 

Figure  22  (Continued) 
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716  FORKAT(T?, 'DISPLACER  NALL ' * T?S, 1 RE 1 2, 4 , T39, I LENGTH  OF  PASSAGE (CM) l 
1#T66# 1PE12,4,T04,6A4,T109, 1 PE  12, 4) 

717  F0RMAT(T2, 'REGENERATOR  SHELL ' # T2S, l PEI 2,4 , T39,  ' EF PICT? VENESS • 

1 #T6S,  IPE12,4,T84,  ' EFFECT  I VENESS  I , T 105,  li»E  12.4) 

710  F0RMAT(T2, 'INSULATION  • . T23, IPE 1 2, 4 ) 

719  F0RmAT(T2#  'REGENERATOR  MATRIX', T29«1PE12,4« 

1 T39# 'LENGTH  OF  FLUI 0 ', T04 ,' LENGTH  OF  FLUID ' /T44 ,« PASSAGE (CM 
1) ' ,T6S,1PE12,4,T89, 'PASSAGE(CM) ' , T 109, 1 PE  1 2,4 ) 

720  FORMAT ( //T21 » 'REGENERATOR  OtSIGN( ' , ?A4, ' ) • /) 

721  FORMAT (//T 1 6, ' MOT  HEAT  EXCHANGER  DFSIGN ( ' , 3A4, ' > ' / ) 

72?  FORMAT (//T1S# ' AMBIENT  MEAT  EXCMANGER  OESIGN( I , 3A4, • ) I /) 

723  FORMAT (T39, 'DYNAMIC  HE  AO ( ATM) ' , T69, l Pt 1 2, 4, T04 , ' OYNAMJC  HEAO(ATM)l 
1,T10S,1PF1?,4) 

724  FORMAT(///T40, 'THERMAL  COMPRESSOR|  RESULTS  OF  THE  ANALYSIS') 

730  FORMAT (///T44, 'THERMAL  COMPRESSOR  DESIGN » //T30, ' GEOMETRY ', T78, 

1 'LOSSES(WATTS) '/T10, 'REOUCEO  VOID  FRAC T ION ' , T43, l PE  1 2 , 4 , TBS, 

2' IDEAL  HEAT  INPUT ' ,T9S, IPE12.4/T10, 'VOLUMETRIC  EFF IC IENC Y « , T4S, 1 PE 
312,4, T63# 'SHUTTLE '#T93, 1PE12.4/T10, 'MASS  FLOP  PER  C YCLE (C) ' , T45, IP 
4E12,4,T65, 'DISPLACFR  CONDUCTION  I , T9S, 1 PE  1 2,4/T 1 0, I MO  END  VOLUME(C 
30',  T4b, 1PE12, 4, T0S, 'CYLINDER  NALL  CONDUCTION' , T93,  IPE12. 4/) 

731  FORMAT(T10, 'AMBIENT  VOIO  VOLUME(CC)',  ' 

• ul  Br  . T4S.IPE12, 4, TBS, 'INSULATION  CONDUCT 

7T0N' ,T9S,  1PE12.4/T10, 'HOT  VOID  VOLUME(CC)',  T45, l PE  1 2, 4, T63,  ' P 
BUMPING' ,T99, 1PE12,4/Tl0, • INTERMEDIATE  VOID  VOLUME(CC)',  T49,lPE 
9j2,4,T69,  ' REGENERATOR  I ' /Ti0, 'CYLINDER  DIAMETER (CM) ' ,T49, 1PE12.4, 

1 T67, ' HF AT  LOSS' ,T09,1PE12,4/)  ’ ' 

732  FORMAT ( T 10, « STROKE (CM) » , T4b, 1 PE  1 2 , 4 , TB7 , 

2 SHELL  CONDUCTION', T99,1PE12,4/T10, 'DISPLACER  ARE A.HOT  END(SQ  CM) I 
3,  T49,1PE12,4,T67, 'MATRIX  CONDUCTION  I , T9S, 1PE 12, 4/T 10, ' 01 SPL ACER  AR 
4M*AMBIENT  fcNO (30  CM)lf 

« . , T49,1PE12,4,T69, MOTAL  HOT  ENO  POWER  INPUT', 

9T99,1PE12,4/T10,'vi  EARANCE  VOLUME. AMBIENT(CC) ' , T49, 1PE12.4/) 

733  FORMAT(T10,  ' 

VOLUME.HOT(CC)',  T49, 1PE12, 4, T76, 'DISPLACER  POWER ( NA 
7TT8) I/T10,  'PERIPHERAL  CLEARANCE  VOLUME(CC)',  T4«,1PE12,4/ 

1T10, 'REGENERATOR  VOID  VOLUME(CC)',  T43, 1PE12.4, T63, 

B'DELTa  P LOSSI,T99,1PE12,4/T10, 'SUCTION  PORT  OPENS(CM) l ,T43, 1PE12. 
94. TBS, 'DELTA  A GAIN',  T9S, 1PE12,4/)  ' 

734  FORMAT  ( T 1 • , 'OISCHARGE  PORT  OPEN 

1S(CM) l ,T4S, 1PE12,4,T69, 'NET  POWER  REQUIREMENT' ,T9S, 1P612. 4) 

END  1 


ROUTINES  CALLEDl 
SORT  , REGDES,  ALOG 

OPTIONS  i/OP | 3 


HXDES 


BLOCK  LENGTH 

MAIN,  9261  (024432)* 

INPUT  12  (000030) 

REG  230  (000714) 

SAVE  400  (001440) 


♦♦compiler 


CORE** 


PHASE  USED  FREE 

DECLARATIVES  01191  08886 
FXECUTABLES  01823  06214 
ASSEMBLY  02487  12190 

Figure  22  (Continued) 
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lUHI 

0002 

0003 

0004 
000b 

0006 

0007 
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000B 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 

0016 


0019 

0020 
0021 


0022 

0023 

0024 


002b 

0026 

0027 

002B 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0036 


C 

C 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


SUBROUTINE  REGDES 

REAL  K4,LR,N,M3A,NRBAR, J0AR,NNU,NTUSAR,NR,KR|,NPR,KRM,KRMM 

C0mM0N/INPUT/TH,TC,TCR,P3,PD,N 

COMMON/ MX / WRB AR, NRBAR 

pnuanii /nrr /uuarr  utu  u in  id  nu  miiriiiai  flDrr.  not  aln  mi 


LU""UN/inruinn,  i fc, 

COMMON/ MX/ WRB AR, NRBAR 

COMMON/RFC/VVRFC,VPH,V3O,lR,OH,OP'i«6U00)#ORrG,ORS,ORM,RRS,TR0, 
lIRfG, IMTX,0#LHC,KRmm 
COMMON/SA VE/WR ( 100) ,P(100) 

DATA  R, CP , NPR, PI /2,079, 5, 193,, 667 ,3,1415926536/ 

REGENERATOR  TMERMAL  design  is  baseo  on  mean  FLOM  parameters, 

J FACTOR  is  FROM  CURVE  FIT  OF  KAYS  AND  LONDON  DATA 


T4«ITH*TC)/A10G(TH/TC) 

U4«3,674F*6«T4**,700 

K4»2,661E"5*T4**,700 

ACbWRFG/LR 

AFRaAC/VFR 

A«4 , ♦ V VMEG/UH 

M3A»,101325*PD*V3D*TCH/(R*TC) 

WRBARb?, *N*M3A 
GBARawRB AR/AC 
NRBAReDH*G8AR/U4 

GO  TO  ( 10,20, 30 ), IMTX 

SPHERES 

10  JBARe,2?l*NRBAR**(«,297) 

HBAReJBARaJ ,30993*CP*GBAR 
GO  TO  40 

SCREENS 

20  JBARa,495*NRBAR**(-,394) 

HBARaJBAR*!, 30993*CP*GB AR 
GO  TO  40 

ROUND  TUBES 

30  Za2,*LR/(NRBAR*NPR*DH) 

IF (Z, IT, ,020)  GO  TO  3b 
Za,43429*AlOG(Z) 

Z*1 ,4307*2*1 ,0 

Z*,7509b72«, 121 1 163* Z*, 02600 129* (2,»Z*Z«1 , ) 
NNU« 1 0 , • * Z 
HBARaNNU*K4/DH 
GO  TO  40 

35  NNUa2,230*Z**(''t391) 

HBARaNNU*K4/DH 

40  NTU0ARe,5*HBAR*A/(WRBAR*CP) 

IF (NRBAR.GT.2000,)  WR1TE<5,750) 

EaNTUBAR/d  .4NTURAR) 

QREGa,S*wRBAR*CP*(TH«TC)*(l,-E) 

Figure  22  (Continued) 
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0039 

0040 


0042 

0043 

0044 
0046 

0046 

0047 
004R 
0049 


0050 

0051 

0052 

0053 

0054 


0055 

0056 

0057 
0056 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 


0067 

0066 


0069 

0070 

0071 

0072 


0073 

0074 


1 


* 


*RITE(5,700)  AC,  APR, A'MRBAR'NRSAR 
WRITF(5,70i)  NTUBAR,E 


REGENERATOR  PRESSURE  DROP  IS  BASED  ON  INSTANTANEOUS  PLOn** 
F FACTOR  16  FROM  CURVF  FIT  OF  KAYS  AND  LONDON  DATA 


60  TO  C50,50,6P),IMTX 


SPHERES  OR  SCREENS 


50  00  55  1*1,100 

C*WR(J)/AC 
NR*DH*ABSCG)/U4 
IF(NR,GT,0,)  60  TO  52 
DPREG ( I ) *0  , 

60  TO  55 

52  Z • ( 2 • * , 4 3429* AL 06 (NR) *5,6990 )/ (»3 ,6990) 

Z**, 055026770*, 561 2S345«Z*, 15975069* (2, «Z«Z*1 ,)♦, 0360671 74* (4 ,*Z*Z 
1 *Z*3,*Z) 

F*10,**Z 

VaR*T4/(P(I )*, 101325) 

DPREGd ) ■ ,5*F* A* V* ASS (6) *6*9,869395*7/ AC 
55  CONTINUE 
60  TO  65 


ROUND  TUBES 


60  DO  61  1*1,100 
G*MR(I)/AC 
NR*DH*ABS(G)/U4 
IF(NR,GT,0.)  GO  TO  62 
DPREG ( I ) *0, 

GO  TO  63 

62  Z*,43429*AL0G(NR«0H/LR) 

Z*(2,*Z. 3, 47721/1, 9246 

Z*1 ,529678* (322540«Z*,0495626l« (2. «Z*Z*1.) 

F*10, **Z/NR 
V*R*T4/Pf I)/, 101325 

63  DPREG ( I ) *,5*F*A*V*A6S(G)«6*9,96939E*7/AC 


CONDUCTION  through  shell 


65  CONTINUE 

60  TO  (60,70), IREG 


regenerator  separate  from  cylinder 


70  DRS*SQRT(4,*AFH/PI) 

MRI TE (5, 702)  DRS 
QRS*KRS*PI*0RS«TRS«(TH*TC)/LR 
60  TO  HP 


regenerator  integral  nith  cylinder 


80  DRSJ*D*2,54 

DRS* SORT (4,*AFR/PI*0RSI*0RSI ) 


Figure  22  (Continued) 
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007b 

0076 

0077 
0070 
0079 
0000 
0001 


0082 

0083 

0084 


0088 

0086 

0087 

0088 

0089 

0090 


0091 

0092 

0093 

0094 


0096 

0096 

0097 

0098 

0099 

0100 


IF(LR-LHC)  100(90(100 
90  QRS>KKS*RI*DRS*TRS*(TH«TC)/LR 
hRITE(5,705)  DR3,D 
CO  TO  110 

100  QR3bKR3*PI*(DR3*DR3I)*TR3*(TH*TC)/LR 
NRI TE (9#  706)  ORSfORSl 
110  CO  TO  (1?0*130,140),ImTX 

MATRIX  CONDUCTION 

SPHERES 

120  RPH«4i*KRMM/(241#*PI) 

ORH«krm*aFR*(TH.TC)/LR 

RETURN 

SCREENS 

130  RHaf25*0H 

ON« 1 0, 4676* (.ALOG10(RH) 1 **1.271) 

RHO0«10#**(#i477#(»4LQglPjpMjj##ig7<j 

KRMaKRMH*HHOW*RHON *00*00*3. 5 A8E«3 
ORM*KRH*AFR*(TH*TC)/LR  j 

RETURN 
ROUND  TUBES 

140  KRM»|9069*KHMM*(VFR.i09Jl) 

ORM*KRH*AFR*(TH«TC)/LR 
600  FORMAT (/SE12,4) 

700  FORMATjTljl^FLON  AREAJSQ  CM)  t,  T45, 1 PFl»,4/Ttt,  l FRONTAL  AREA  (80  CM) 
1T45,1PE12,4)C  ’ * T46,1PE12,4/T10,»MFAN  REYNOLDS  NUMBER > , 

t"‘n5fer  me, 

702  FORMAT (Tit* l DIAMETER (CM) I , T45. tPEiz  4) 

0SiE°rM5ii^?^?ci?f?54s:r^iS.i27ETcR<£H,,#T45'ipci2-4/Tj0*,i^i'>e  sm 

750  RETURNt//TI0* '**********  FL°*  IS  TURBULENT  **********1//) 

end 


SFORTRN  BIlMXDES 


ROUTINES  CALLEDi 
4L0G  , A63  , 3URT 

OPT ION3  a/0P|3 


ALOG10 


BLOCK 

REGDES 

INPUT 

HX 

REG 

SAVE 


LENGTH 

1444  (0e6910)* 

12  (000030) 

A (000010) 

230  (000714) 

400  (001440) 

Figure  22  (Continued) 
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U00 1 
0002 
1003 
( 004 
(005 

( 00b 
( 007 
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SKrDs:i?5.sass;wf;£!«;stfr£%rt75»D,^“,L"L'r'B,,^'*'B^^* 

DIMENSION  OPHX(100),OPFP(100) 

COMMON/HX/PRBAR, NRBAR 
COMMON/INPUT/TH, TC»  TCR,PS,PD,N 
COMMON/S A VE/WR( 100 ),P()00) 

OAT A «»CP,NPR,PI/2,0P0,0, J93, i66T#3j l4 jBg269j6/ 

THERMAE  DESIGN  OF  THE  HEAT  EXCHANGER  IS  BASED  ON  MEAN  FLOP  PROPEPT 

U*3,674E-6*T**,700 

TWG0«2,*0INBAR*(1,.E)/E/(pRBAR*CP) 

RHO", J0|325*PS/R/T 

AFPiWRRa'/SQRT (OTNH*RHO/9tBb939E»7 ) 

vfp*afp«lfp 

VHX«VV»VCL"VFP 

NTU«AlOGU,/(l,.E)J 

WRITE (5, 700)  TPGO,RHO,AFP,VFP,VHX,NTU 

IF(VHX,GT,0,)  GO  TO  20 
WRITE(5,501) 

RETURN 

20  BAN«(0*2,54)*PI 

CALCULATE  THE  PRESSURE  DROP  IN  THE  FLUID  PASSAGE 


SQAiSQRT(AFP) 

DO  25  1*1,100 
AmR*ABS ( wR ( I ) ) 

I F ( A RR , Nt , 0 , ) GO  TO  22 
DPFP(I)*0, 

GO  TO  25 

22  NR*AMR/(2(*u«3QA) 

Fa  1 ( 025 3* NR * • (• ,9422 ) 

__  I>g4«**,86039E*7*F*(LFP/3QA/RHO)*PR(I)*APR/AFP/AFP 

ca  CONTINUE 

PRESSURE  DROP  IN  THE  HEAT  EXCHANGER  IS  BASEO  ON  INSTANTANEOU*  FLOP 

GO  TO  (30,40) , IHX 

ANNULAR 

30  NRB ARa2 ( aPRBAR/BAN/U 

IF ( NRBAR , GT ,2000, ) PRITE(5,750) 

N8T*7(54/NRBAR/NPR 

L*SURT(,5*NTU*VHX/NST/BAN) 

CaSQRT(?,*VHX*N8T/BAN/NTU) 

0H*2, *C 

PRI TE (5, 701  ) NRBAR, NST, BAN, C,L,E 
DO  35  1*1,100 
G*PR(I)«L/VHX 
V*R*T/P(I)/, 101325 

35  DPHX(I)*48,0*L*U*V*G*0,B6939E*7/OH/OH 
Figure  22  (Continued) 
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0042 


RETURN 


C MULTIPASSAGE. 

C 

0043  40  BPeSQRT(4f*2,B9*U*VHX/(NT.U*NPR*MRBAR)) 

0044  NP«VMX/CHP*BP*LP) 

0045  NRBAReWRBAR/ (NP*BP*U) 

0046  IMNRBAR,GT,2000, ) mRITE(5,750) 

0047  DHeBP 

0048  NRI TE (5,  702 ) BP,NP,LR 

0049  00  45  I ® l , 100 

0050  G»mH(I)/CNP*BP*8P) 

0051  VeR*T/P(I)/,10l3?5 

0052  4b  DPHX(I)e28,4*L*U*V*G*9,86939E-7/0H/0H 

0053  RETURN 

0054  500  F0RMATC/2X, IHX' ,0CIPE12,4))  „ ...  _ .... 

0055  501  FORMAT (// ' ***  HEAT  EXCHANGER  VOLUME  IS  ZERO  OR  NFGATlvE  ***'//> 

0056  700  FORMAT  C T 1 0*  'EXIT  DELTA  T (K) ' , T45, 1PE12,«/T10,  GAS  DEMSITY(G/C 

1C)',  T45, 1PE12.4/T10, 'FLUID  PASSAGE  AREA (SQ  CM)', T 45,1 PE  12, 4/ 

1TI0,  'FLUID  PASSAGE  VOLUME(CC)',  T45, 1PF 1 2 , 4/Ti 0, ' HE AT  £*j*HANGE 

1R  VOLUME (CCI',T45, IPE12,4/T1«, 'HEAT  TRANSFER  UNITS '#J45,i PF12, 4) 

0057  701  FORMAT(T10, 'MEAN  REYNOLDS  NUMBER ', T45, 1PE12,4/T10, ' MEAN  STANTON  NU 

lMBER' ,T45,  iPFl2,4/T10, 'PERIMETER  OF  ANNuLUS (CM) ' , T45, iPE 1 2, 4/ 

1 T 1 0 , ' ENTRANCE  NIDTH(CM)',T45,1PF.12,4/T10,»LENGTH(CM)'*T45,  IPE  12,4/ 
1T10, 'EFFECTIVENESS' ,T45,1PE12, 4)  . . .. 

0056  702  FORMAT(T10,  IPASSAGF  N10TH(CM ) » , T45, 1 PE  1 2, 4/T 1 0 , ' NUMBER  OF  PASSAGES 

1 ' ,T45, 1PF12,4/T10, 'LENGTH  OF  PASSAGE (CM) ', T45, IPE 12, 4) 

0059  750  FORMAT(//TI0, ***********  FLOW  IS  TURBULENT  *********♦'//) 

0060  FND 


ROUTINES  CALLFOl 
SORT  , ALOG  , ABS 

OPTIONS  «/0Pl3 

block  LENGTH 

HXDES  1112  (P04260)* 

HX  4 (000010) 

INPUT  12  (000030) 

SAVE  400  (001440) 

• •COMPILER  CORE** 

P4ASf  USED  FREE 

DECLARATIVES  00622  09415 
EXECUTABLES  0129/  08740 
ASSEMBLY  01659  13018 


Figure  22  (Continued) 
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This  report  presents  the  results  of  a design  study  program  on  a rotary 
reciprocating  compressor  which  uses  heat  as  the  source  of  driving  energy.  The 
objective  of  the  study  program  was  to  determine  the  feasibility  of  using  a ther- 
mal compressor  In  a spaceborne,  Cryogenic  refrigeration  system.  The  report  ex- 
plains the  principles  of  operation  of  a thermal  compressor,  and  presents  the 
equations  which  describe  Its  behavior.  It  also  describes  a computer  program  con- 
structed to  manipulate  these  equations.  The  results  of  a number  of  runs  made 
with  this  program  are  discussed,  and  the  design  of  a thermal  compressor  to  meet 
a specific  set  of  performance  requirements  Is  described.  Finally,  a comparison 
Is  made  between  two  types  of  spaceborne,  closed-cycle,  cryogenic  refrigeration 
systems — one  Incorporating  a thermal  compressor  and  the  other  Incorporating  an 
electrically  driven  compressor. 

It  Is  concluded  that  a thermal  compressor  Is  a viable  alternative  to  an 
electrically  driven  compressor  In  refrigeration  systems  where  heat  la  preferred 
over  electricity  as  the  energy  source.  For  the  specific  case  studied,  the  re- 
frigeration system  Incorporating  a thermal  compressor  required  92  watte  of 
electrical  power  and  1,064  watts  of  heat  at  1,250#F;  while  the  system  Incorpora- 
ting an  electric  compressor  required  587  watts  of  electrical  power.  The  weights 
and  sizes  of  the  two  refrigeration  systems  were  comparable.  No  serious  techni- 
cal problems  are  foreseen  In  building  a thermal  compressor . 
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